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ABSTRACT 

We present the results of an ultradeep, narrowband imaging survey for Lyman-continuum (LyC) 
emission at z ~ 3 in the SSA22a field. We employ a custom narrowband filter centered at A = 3640 A 
(NB3640), which probes the LyC region for galaxies at z > 3.06. We also analyze new and archival 
NB4980 imaging tuned to the wavelength of the Lya emission line at z = 3.09, and archival broadband 
B, V, and R images of the non-ionizing UV continuum. Our NB3640 images contain 26 z > 3.06 
Lyman Break Galaxies (LBGs) as well as a set of 130 Lya emitters (LAEs), identified by their excess 
NB4980 flux relative to the BV continuum. Six LBGs and 28 LAEs are detected in the NB3640 
image. LBGs appear to span a range of NB3640— R colors, while LAEs appear bimodal in their 
NB3640— R properties. We estimate average UV to LyC flux density ratios, corrected for foreground 
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contamination and intergalactic medium absorption, finding (Fuv I FLyc)cor? = ii -"-5 4 
implies a LBG LyC escape fraction f^jf ~ 0.1, and (F uv /F Ly c)^f = 2 - 2 ~oi- The strikingly blue 
LAE flux density ratios defy interpretation in terms of standard stellar population models. Assuming 
{Fuv/FL V c)carr applies down to L = 0.1L*, we estimate a galaxy contribution to the intergalactic 
hydrogen ionization rate that is consistent with independent estimates based on the Lya forest opacity 
at z ~ 3. If we assume that (Fjjv I FLyc)corr holds at the faintest luminosities, the galaxy contribution 
significantly exceeds that inferred from the Lya forest. We interpret our results in terms of a model 
where LyC photons escape over only ~ 10 — 20% of solid angle. When advantageously oriented, a 
galaxy will exhibit a low UV to LyC ratio, an effect enhanced for more compact galaxies. This model, 
however, does not adequately explain the extremely blue NB3640— R colors measured for some LAEs 
in our sample. Further follow-up study of these faint LAEs is crucial, given the potentially important 
contribution similar objects make to the process of reionization. 

Subject headings: galaxies: high-redshift - intergalactic medium - cosmology: observations - diffuse 
radiation 



1. INTRODUCTION 

The thermal state of the intergalactic medium (IGM) is 
governed by its spatial density distribution, the spectral 
shape and intensity of the metagalactic ionizing back- 
ground, and the distribution of the sources of the back- 
ground. Identifying these sources constitutes a criti- 
cal step towards understanding the physics and ther- 
mal history of the IGM. QSOs are efficient producers 
of far-UV radiation and dominate the Lyman-continuum 
(LyC) emissivity that maintains the ionization state of 
the I GM up to the pea k of AGN activity at redshift 
z ~ 2 (|Cowie et al.ll2009|) . Beyond z ~ 2, however, their 
contribution decr eases drastically wit h their decreasing 
number density (jHopkins et al.l (20071 ) . As it has been 
demonstrated that reionization was completed by z > 6 

1 Based, in part, on data obtained at the W.M. Keck Obser- 
vatory, which is operated as a scientific partnership among the 
California Institute of Technology, the University of California, 
and NASA, and was made possible by the generous financial sup- 
port of the W.M. Keck Foundation. 

2 Based, in part, on data collected at the Subaru Telescope 
and obtained from the SMOKA archive, which is operated by 
the Astronomy Data Center, National Astronomical Observa- 
tory of Japan. 

3 David and Lucile Packard Fellow 



(|Fan et al.l [200llBecker et all 120071) . it is generally as- 
sumed that the ionizing flux required for reionization was 
produced in star-forming galaxies. Direct measurements 
of ionizing radiation escaping from galaxies at z > 6 are 
not possible because of the large optical depth to LyC 
photons in the high-redshift IGM. However, the measure- 
ment can in principle be made at lower redshift where the 
IGM is less opaque. 

Although observations of star-forming galaxies at rel- 
atively low redshift (0.5 < z < 2) have failed to de 



tect escaping LyC photons (e.g.. iMalkan et al.l 
120071: iGrimes et all 120091: ISiana et ail 



Siana et all 



2003 



2010; 



Bridge et al(i20Tq) , there has been much recent progress 



in the search for LyC emission at higher (z ~ 3-4) red- 
shift. Building on the initial detection by ISteidel et all 
(120011) in a composit e spectrum of 29 LBGs at (z) = 3.4, 
Shapl ev et alJ (|2006[ ) presented the first spectroscopic de- 
tections of LyC emission from individual objects. In 
a sample of 14 objects with deep Keck/LRIS spectra, 
two exhibited significant LyC emission, implying f^ s v c > 
0.15. Deep Keck/LRIS spectroscopy for a z ~ 3 sam- 
ple larger by an order of magnitude is forthcoming, and 
can be used to generalize the initial spectroscopic results 
(Steidel et al, in preparation). 
An alternative technique for studying the escape of 
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LyC emission consists of narrowband imaging. Measur- 
ing the density of LyC radiation with imaging requires 
a filter that only transmits flux below the Lyman limit, 
and also one that is relatively narrow. A filter cover- 
ing too broad a range below the Lyman limit will indi- 
cate more about the statistics of IGM absorption than 
the absorption intrinsic to galaxies. Specifically, a fil- 
ter width comparable to or larger than than the wave- 
length range corresponding to one LyC photon mean free 
path will typically contain significant absorption due to 
intervening Lyman limit systems, thus making it more 
difficult to infer the amount of escaping ionizing radia- 
tion in the immediate vicinity of a high-redshift galaxy. 
Current estimates of the mean free path at z ~ 3 are ~ 
70 — 85 Mp c (proper) , assuming the current concordance 
cosmology (iFaucher-Giguere et al.l2008l : IProchaska et al.1 
12003 ISongaila fc Cowiell2010h . This range in mean free 
path corresponds to a rest-frame wavelength interval of 
~ 850 — 910 A. Since standard broadband U filters probe 
a larger wavelength baseline than this, it is advantageous 
to adopt a custom setup, including a narrowband filter 
whose transmission curve probes within the rest wave- 
length range corresponding to one mean free path. 

Key for maximizing the efficiency of narrowband imag- 
ing in a fixed rest-frame wavelength range is the existence 
of a large sample of galaxies in a single field at roughly the 
same redshift. In this respect, the SSA22a field serves as 
an ideal target for narrowband imaging searches for LyC 
emission at z ~ 3. This field has been shown to contain a 
significant overdensi ty of Lyman Break Galaxies (LBGs) 
at z = 3.09 ± 0.03 (|Steidel et all 119981 ). In addition to 
the 27 LBGs spectroscopically confirmed to lie within the 
z = 3.09 redshift spike, narrowband imaging tuned to the 
wavelength of Lya emission at z = 3.09 has been used 
to identify large samples of Lya emitters (L AEs) and ex- 
tended Lyct "blobs" at the same redshift (ISteidel et al.l 
12000b iMatsuda et alJl20"ol iHayashino et al II2004T ). 

Escaping LyC emission at z > 3 detected by n arrow - 
band imaging was first reported by llwata et a l. (2009). 
Using Subaru/Suprime-Cam and a "NB359" filter, with a 
central wavelength 3590A, and FWHM of 150A, these au- 
thors searched for LyC emission for sources at z > 3.06, 
and reported detections for 7 LBGs and 10 LAEs. In 
some cases there are spatial offsets of several kpc between 
the centroids of ionizing and non-ionizing UV emission, 
which may distinguish amon g different models for the 
escape of LyC emission (e.g.. iGnedin et al.ll200ll . Fur- 
thermore, in some of the fainter LAEs, the high apparent 
ratio of escaping ionizing to non-ionizing UV radiation 
challenges standard models for the intrinsic spectral en- 
ergy distribution of star-forming galaxies (jlnoud 120101 : 
llnoue et alj l20lif ). In a parallel and independent ex- 
periment, we have used Keck/LRIS and a narrowband 
"NB3640" filter, similar to the Suprime-Cam NB359 fil- 
ter, in order to collect even deeper images of the LyC 
region for z ~ 3 galaxies in the SSA22a field. We have 
also obtained significantly deeper imaging tuned to the 
wavelength of Lya, enabling more robust identification 
of LAEs, down to fainter flux limits. 

In this paper, we present the results of our narrowband 
imaging survey for LyC emission at z ~ 3 in the SSA22a 
field. 321 describes the observations and data reduction. 
In §3, we present narrow and broadband photometric 



measurements, the identification of a sample of LAEs 
based on our new, deeper data set, and the method used 
for matching apparent LyC detections with known z ~ 3 
targets. We also discuss potential emission contamina- 
tion along the line of sight in our narrowband imaging. 
§4 contains our results, including estimates of the true 
ratios of ionizing to non-ionizing flux densities for our 
targets, the corresponding escape fractions, and the im- 
plied space density of ionizing radiation. In §5, we dis- 
cuss the cosmological implications of these results, and 
offer a concluding summary in §6. Throughout the pa- 
per we employ the AB magnitude system, quote flux and 
luminosity densities in units of F v and L v , respectively, 
and assume a cosmology with f2 m = 0.3, fl\ = 0.7, and 
H = 70 km s- 1 Mpc" 1 . 

2. OBSERVATIONS AND DATA REDUCTION 
A custom narrowband filter was used for LyC obser- 
vations, with central wavelength 3635 A and FWHM 
100 A, hereafter referred to as NB3640. The filter was 
originally designed to detect Lya emission at z = 2 
(Andrew Blain, private communication), but its spec- 
ifications are such that no significant contamination 
(< 1%) longwards of the Lyman limit (912A) will en- 
ter the filter for redshifts z > 3.06. Furthermore, the 
narrow NB3640 bandpass probes the wavelength region 
well within one mean free path below the Lyman limit 
at t he redshift for the majority of our target galax- 
ies (iFaucher-Giguere et"afl 120081: iProchaska et all 120091 : 
ISongaila fc Cowiel 120101 ) thus minimizing the correc- 
tion for intergalactic medium absorption (IGM ) that 
is required to infer f^jf (jShaplev et al.l 120061 ) . Fig- 
ure [1] shows the transmission curve of the NB3640 fil- 
ter o verlaid on a compos ite z = 3.09 LBG spectrum 
from IShapley et al.1 (|2006[ ) . Also shown in Figure [1] is 
the transmission curve of an additional custom narrow- 
band filter tuned to the wavelength of Lya at z = 3.09, 
the SSA22a spike redshift. This filter has central wave- 
length 4985 A and FWHM 80 A, hereafter referred to as 
NB4980. 

We obtained NB3640 imaging data using the blue 
side of the Low Re solution Imaging Spectrometer (LRIS; 
ISteidel et af1l2004f ). The NB3640 imaging data were col- 
lected on 6-7 August 2005 and 2-4 June 2008 (UT), for 
a total of 13.2 hours of integration. During the August 
2005 run, the red side of LRIS was simultaneously used 
to collect NB4980 data, for a total integration time of 
9.4 hours. The blue side (LRIS-B) detector has a plate 
scale of 0'.'135, while the red side (LRIS-R) has 0'.'211 
pixels. The "d460" dichroic beam splitter was used to 
send light bluer than ~ 4600 A towards the LRIS-B, 
while longer wavelengths were directed towards LRIS-R. 
Conditions were photometric during both observing runs, 
with seeing in the final NB3640 and NB4980 stacked im- 
ages of 0'.'80. The pointing center of the LRIS imaging 
field was a =22:17:28.8, S =+00:14:36, and the area of 
the final trimmed stacks is 5.5' x 7.6'. Contained within 
this footprint are 26 of the 33 LBGs in SSA22a with 
spectroscopic redshifts z > 3.06, such that the NB3640 
probes bluewards of the Lyman limit. There are also 
two z > 3.06 QSOs, and 41 LBG photometric candidates 
without spectroscopic redshifts. The images were taken 
at a sky position angle of 9 = 180°, with the long axis 
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Fig. 1. — Narrowband filter transmission curves overlaid on an LBG composite spectrum from [Shaplev et'a l. (2006), with 
(z) = 3.09. The NB3640 filter probes the LyC region for galaxies at z > 3.06, while the NB4980 filter contains the Lya feature 
for galaxies at z ~ 3.09, the redshift of the significant galaxy overdensity in the SSA22a field. Dotted lines mark both the 
Lyman limit and Lya feature at z = 3.09. 



of the image oriented along the North-South direction. 
The LRIS-B 4K x 4K CCD mosaic detector consists of 
two 2K x 4K chips with slightly different (~ 20%) quan- 
tum efficiency at 3640A, and the pointing center and sky 
position angle were chosen to maximize the number of 
z > 3.06 targets falling on the blue side chip with better 
quantum efficiency. 

Individual narrowband exposures were reduced using 
standard procedures of flatfielding, background subtrac- 
tion, and removal of cosmic rays. For both NB3640 and 
NB4980 imaging, flatfields were constructed from de- 
registered images of the twilight sky. For imaging in each 
narrowband filter, tasks from the IRAF "images" and 
"mscred" packages were used to solve for the world co- 
ordinate system of each exposure, sample each exposure 
onto a reference world coordinate system, and shift and 
combine all individual exposures into a final stacked im- 
age. During combination of individual exposures, which 
ranged in airmass from 1.06 — 1.37, scaling was applied so 
that each exposure effectively was observed at the mini- 
mum airmass. Furthermore, LRIS-B images were resam- 
pled onto the 0'.'211 LRIS-R pixel scale. Narrowband 
images were calibrated onto the AB system using obser- 
vations iof_sj2e£troj2hotometric standard stars from the 
list of iMassev etaLl (jl988T ). NB3640 and NB4980 AB 
magnitudes were also corrected for a Galactic extinction 
of E(B - V) = 0.08, based on I RAS 100/zm cirrus emis- 
sion maps (ISchleeel et al.lH998f) . 

We augmented the LRIS data set with both narrow- 
band and bro adband archival imag ing data taken with 
Suprime-Cam (Miyazaki ct al. 2002) on the Subaru Tele- 
scope. Additional data include narrowband images taken 
through a filter with central wavelength 4977 A and 
FWHM 80 A, i.e. almost identical to the one we used 
for Keck/LRIS observat io ns, and describe d prev iously 
in IMatsuda et all (I2004D. lHavashino et al.l (|2004D . and 
llwata et all (|2009D . These "NB497" images were ob- 
tained on 8-9 September 2002 (UT), for a total of 7.2 
hours of integration time, and an image quality in the 
stack of 0'.'76 FWHM. Broadband B- and U-band im- 



ages were obtained on 21 September 2003 (UT) (PI: Hu) 
with integration times of 1.1 and 0.8 hours, respectively, 
and respective image quality of 0'.'76 and 0'.'60. Finally, 
we analyze i?-ba nd imaging from 20 October 2001 (UT) 
(|Hu et al.ll2004D . with an integration time of 0.8 hours 
and seeing of 0'.'56. All Suprime-Cam da ta were re- 
trieve d from the Subaru archive, SMOKA (Ba ba et al.1 
120021), and reduced with the SDFRED reduction pack- 
age (|Yagi et ^12001 lOuchi et al.ll20"ol . Subaru B, V- 
band, and R images were flux calibrated u sing a combina- 
tion of spectrophotometric and Landolt (jLandoltl 119921) 
standard stars. A summary of both Keck/LRIS and 
Subaru/Suprime-Cam observations is presented in Ta- 
ble [Q 

Given the similarity of the Lya narrowband filters used 
for LRIS and Suprime-Cam imaging, in terms of cen- 
tral wavelength and FWHM, we combined NB4980 and 
NB497 images, scaling them to a common zeropoint, 
and weighting them according to their respective depths. 
This combined NB4980 image is ~ 0. 4 mag deeper 
than t he Subaru-only image presented in Matsud a et al.1 
(2004) . The selection of Lya emitters (LAEs) requires a 
continuum image probing the same effective wavelength 
as the narrowband. Since the combination of B and V 
filters straddles the NB4980 bandpass, a linear combi- 
nation of the two filters suffices for the matched con- 
tinuum image. We constructed a BV continuum im- 
age by first smoothing the U-band image so that the B 
and V image PSFs were matched, and then scaling the 
PSF-matched images to the same zeropoint. We then 
executed: BV = 0.44B + 0.56U, with the weights de- 
termined by the relative wavelengths of the broadband 
and NB4980 filters. We note that the linear combination 
of B and V images d escribed here differs from that in 
IMatsuda etldl (|2004D (i.e., {2B + U)/3), and better re- 
flects the relative wavelengths probed by the B, V, and 
NB4980 filters. Finally, we smoothed the BV continuum 
image so that the PSF matched that of the AB4980 com- 
bined image. As described in §3.2[ selection of LAEs is 
based on the NB4980 magnitudes and BV- AB4980 col- 
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TABLE 1 
Description of Observations 



Filter Telescope/Instrument Seeing FWHM (") Exposure (s) 

NB3640 Keck/LRIS-B 0.80 47636 

NB4980 Keck/LRIS-R 0.80 33880 

Subaru/Suprime-Cam 0.76 25800 



B Subaru/Suprime-Cam 0.76 3927 

V Subaru/Suprime-Cam 0.60 3000 

R Subaru/Suprime-Cam 0.56 3000 



ors of objects. However, for examining the spatial distri- 
bution of Lya emission (and absorption), isolated from 
the local continuum, we also constructed a continuum- 
subtracted Lya image, hereafter referred to as LyA. 

To estimate the noise properties and depth of the im- 
ages from which photometry was measured (NB3640, 
LRIS+Suprime-Cam NB4980, W-continuum, and R), 
we evaluated the counts in apertures distributed ran- 
domly in regions free of detected s ources. To facilitate 
comparison with llwata et al.l (|2009D . we used 1'.'2 diam- 
eter circular apertures. The standard deviation of the 
counts in empty apertures yields the associated la mag- 
nitude limit, from which an estimate of the 3er limiting 
magnitude can be derived. We find 3cr limiting mag- 
nitudes of 27.90, 27.90, 28.40, and 27.60, respectively, 
for the NB3640, combined LRIS+Suprime-Cam NB4980, 
BV^-continuum, and R images. Our NB3640 image is 
therefore roughly 6 deeper th a n the analogous NB359 
image presented in llwata et al.l ([20091 ) . As described in 
H3.ll this simple estimate of the noise properties of our 
images does not take into account the full systematic 
uncertainty associated with source detection and eval- 
uation. Monte Carlo simulations replicating our detec- 
tion and measurement process with artificial galaxies of 
known magnitude indicate larger photometric errors than 
inferred from the stated depths above. All error bars 
presented in later sections are based on the more conser- 
vative Monte Carlo simulations. 

Finally, for detailed rest-frame UV (non-ionizing) con- 
tinuum morphologies of our SSA22a targets, we assem- 
bled archival Hubble Space Telescope (HST) Advanced 
Camera for Surveys (ACS) imaging of the field, taken 
with the F814W filter. These data i nclude 3 point- 
ings f rom the Gemini Deep Deep Survey ([Abraham et al.l 
120071 P.I.D. 9760), with 5-orbit exposure times, and 3 
pointi ngs from the program described in IGeach et al.l 
(2009) (P.I.D. 10405), with 3-orbit exposure times. Of 
the 34 targets with NB3640 detections described in H3AI 
25 have deep ACS coverage. 

3. ANALYSIS 

3.1. Photometric measurements and uncertainties 

Source identification and ph otometry were performe d 
using SExtractor version 2.5.0 (jBertin fc Arnoutslll996[) . 
We first smoothed the NB3640 image using a circular 
Gaussian kernel with a = 1 pixel = 0'.'211. We then ran 
SExtractor in "dual-image mode" , using the smoothed 
image to detect sources though the measurements were 
performed on the unsmoothed image. We found this 
method produced superior results in terms of detecting 



faint sources while not producing overly-large apertures, 
compared to the use of SExtractor's built-in filtering op- 
tion. Similarly, the NB4980 image was used to identify 
sources and define isophotes for measurements in the BV 
image for the purpose of identifying LAE candidates (see 
below) . 

Due to the details of the radiative transfer and sources 
of the corresponding photons, any detected Lya emission 
and escaping LyC flux will not necessarily be co-spatial 
with either each other or with the bulk of the rest-frame 
UV flux in a given galaxy. Therefore, single-image mode 
was employed for measurements in the NB4980 and R 
images with no prior used to define apertures, and we 
report "Kron-like" (i.e., SExtractor MAG.AUTO) total 
magnitudes for these and the NB3640 images. All re- 
ported colors are determined by taking the differences 
of these magnitudes. This method allows more flexi- 
bilty (e.g., in avoiding contamination from neighboring 
sources) and decreases the photometric uncertainty, as 
compared to using apertures of fixed shape and size large 
enough to capture all of the flux from the most extended 
sources. Furthermore, we restrict "detections" to magni- 
tude ranges for which any systemics in the photometry 
are small compared to the random error (see below). 

In addition to the estimates of the depth reached in 
each image described above (@, we ran simulations of 
our method in order to obtain a detailed description of 
the combined systematic and statistical uncertainties in 
our measurements. We placed one hundred simulated 
sources of known magnitude into the NB3640 image, the 
R image, and of known magnitude and color simultane- 
ously into the NB4980 and BV images. The simulated 
sources were placed at random positions chosen to avoid 
both real and previously placed simulated sources. SEx- 
tractor was then run on these simulated source-added im- 
ages in the same manner as above (i.e., in single- image 
mode on R and NB4980 and dual-image mode for de- 
tection/measurement on smoothed/unsmoothed NB3640 
and NB4980/£?y) and magnitudes were recorded for all 
recovered simulated sources. This process was repeated 
until a large enough total number (typically ~ 50000) of 
simulated sources had been recovered to populate each 
bin in magnitude (and color) to a statistically significant 
degree. The average differences and standard deviations 
of the recovered from input magnitudes then define any 
systematics and the uncertainty as a function of magni- 
tude (and color). We list the systematics and uncertain- 
ties in NB3640 and R in Tabled and present those for 
NB4980 and W-NB4980 in Tabled 
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TABLE 2 

Uncertainties in simulated photometry: R and NB3640. 



magnitude range a AR b cr R c ANB3640 b cr N B3640 c 



23.0 - 


23.5 


0.00 


0.05 


-0.01 


0.03 


23.5 - 


24.0 


0.00 


0.08 


-0.01 


0.04 


24.0 - 


24.5 


-0.01 


0.11 


-0.02 


0.06 


24.5 - 


25.0 


-0.01 


0.16 


-0.03 


0.10 


25.0 - 


25.5 


-0.02 


0.23 


-0.04 


0.15 


25.5 - 


26.0 


-0.01 


0.29 


-0.06 


0.20 


26.0 - 


26.5 


0.01 


0.37 


-0.06 


0.29 


26.5 - 


27.0 


0.13 


0.40 


-0.03 


0.38 


27.0 - 


27.5 


0.43 


0.55 


0.03 


0.49 



a Recovered range of R or NB3640 magnitude for simulated galax- 
ies. 

b Average value of the difference between measured and input mag- 
nitudes. Significant departures from zero imply systematic biases 
in the photometry. 

c The standard deviation in the difference between measured and 
input magnitudes. 

3.2. Lya emitter candidates 

The wavelength of redshifted Lya falls within our 
NB4980 filter for sources at 3.054 < z < 3.1200 Fig- 
ure [2] shows isophotal BV— NB4980 color as a function 
of total NB4980 magnitude for all sources detected in 
the NB4980 image. As the isophotes defined in the 
NB4980 image were used for the BV photometry, source- 
matching was not an issue. Note that obje cts identified 
by S Extractor within Lya "blobs" (e.g., iSteidel et al.l 
2000) were removed from the catalog. The solid line 
indicates our formal lower-limit in BV— NB4980 color. 
For sources bright in NB4980, the limit in color depends 
primarily on the detection limit BV = 27.5. At fainter 
narrowband magnitudes, the color limit is determined 
empirically from the Monte Carlo simulations described 
above. 

Following lSteidel et al.l (|2000j ) , we identify as LAE can- 
didates sources having W-NB4980 > 0.7. The dotted 
lines in Figure [2] indicate the ±3<r scatter of the main dis- 
tribution of sources, which is dominated by objects out- 
side of the redshift range 3.054 < z < 3.120. The right- 
hand axis of Figure [2] shows the Lya rest-frame equiva- 
lent width, EW , implied by the W-NB4980 color for 
sources in this redshift range. The limit BV— NB4980 

> 0.7 (corresponding to an observed Lya equivalent 
width ~ 80A) is above this scatter down to NB4980 ~ 26. 
A total of 110 sources meet the criteria BF-NB4980 

> 0.7 and NB4980 < 26, which we use to define our main 
LAE candidate sample. For reference, our simulations in- 
dicate the uncertainty in color is a W-NB4980 — 0.2 for 
W-NB4980 = 0.7 and NB4980 = 26. We also consider 
sources with 26 <NB4980 < 26.5 with the more conser- 
vative color-cut of BV— NB4980 > 1.2 as an additional 
"faint" LAE sample. Twenty sources meet these crite- 
ria. As these sources arc selected using different criteria, 
they are excluded from any following calculations using 
the ensemble of LAEs, unless otherwise noted. The limits 
defining these sample are shown in Figure [5] with dashed 
lines. 

To test the robustness of the LAE selection tech- 

4 Our NB3640 filter is opaque to non-ionizing UV flux above 
z > 3.06 and transparent to only negligible levels of non-ionizing 
flux at 2 ~ 3.054. 



nique, we compared the equivalent widths implied by 
the W-NB4980 colors for the 18 LBGs in the field 
with 3.05 < z < 3.12 (excluding three within Lya 
"blobs" ) to those determin ed from existing LRIS spectra 
(jShaplev et al.1 12003L 120061 ). The equivalent width val- 
ues determined from the two methods correlate strongly. 
The photometrically determined widths are on average 
larger, which may be due to slit-losses in the spectra 
as the spatial extent of the Lya emis sion tends to be 
larger than that of the r est-frame UV (jHayashino et al.l 
120041 : ISteidel et al.ll20ilT ) . Importantly, the spectra of all 
five LBGs selected as LAEs exhibit Lya in emission, and 
no LBGs with spectroscopic observed equivalent widths 
greater than 80A are missed. Additionally, we verified 
that none of the 45 galaxies with spectroscopic redshifts 
z < 3.05 in the field were selected as LAEs. We also ran 
SExtractor on the LyA (NB4980- W difference) image; 
all of the 130 LAE candidates were recovered as detec- 
tions. 

3.3. Object matching 

For each LBG, we compiled the positions and photo- 
metric measurements for any SExtractor detection in the 
NB3640 image within a 2" (corresponding to 15 kpc at 
z = 3.09) offset from the center of the i?-band aperture. 
For each LAE candidate (hereafter, LAE), we retained 
any detections in the NB3640 and R image within 2" of 
the NB4980 aperture center. We visually inspected each 
matched source in all available bands, including the deep 
BV and the high image quality _R, and removed any obvi- 
ous false matches corresponding to neighboring sources. 
In rare cases, the dcblending produced by SExtractor 
was unsatisfactory and the photometry was recomputed 
using the IRAF task phot. 

Six of the 26 LBGs were detected in our NB3640 image. 
We retain the NB3640 detection near C16 despite the 
large (l'/9) offset, as there are no neighbors detected in 
any of the other bands coincident with the position of 
the NB3640 flux. Of the 110 main sample LAEs, 27 
were detected in NB3640 and 72 were detected in R with 
R < 27, which we adopt as our formal detection limit 
based on our simulations discussed in §3.11 As with C16, 
a small number of detections near LAEs having fairly 
large (~ 1'.'5 - 2") offsets are retained as the LAE is 
the closest source detected in any of the other bands. 
We account for the possibility of mismatched NB3640 
detections in £ 13.51 One of the 20 faint sample LAEs was 
also detected in the NB3640 image, and two detected in 
the i?-band image with R < 27. 

3.4. LBGs and LAEs with NB3640 detection 

Table |4] lists the coordinates, redshifts based on Lya 
emission and interstellar absorption (when available), R 
and NB3640 magnitudes, and spatial offsets of the R 
and NB3640 emission centroids (Ar) for the six LBGs 
with NB3640 detect ions. Of the two LBGs reported by 
iShaplev et al.l (|2006f ) to exhibit LyC flux in their spec- 
tra, we detect one (C49), consistent with the results of 
llwata et al.l (|2009h . The 28 (27 main and one faint sam- 
ple) LAEs with NB3640 detections are marked with cir- 
cles in Figure [5J and the coordinates, photometry, Lya 
EW , and A R and LyA-NB3640 (A LyA ) spatial offsets 
are presented in Table [5] We also note that five of the 
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TABLE 3 

Uncertainties in simulated photometry: NB4980 and _BV-NB4980. 



NB4980 ; 



'-NB4980 a 


ANB4980 b 


CNB4980 C 


A(W-NB4980) b 


"BV-NB3 


0.0 - 0.2 


0.01 


0.03 


0.00 


0.03 


0.2 - 0.4 


0.01 


0.04 


0.00 


0.03 


0.4 - 0.6 


0.01 


0.04 


0.00 


0.03 


0.6 - 0.8 


0.01 


0.03 


-0.00 


0.03 


0.8 - 1.0 


0.01 


0.04 


-0.00 


0.04 


1.0 - 1.2 


0.01 


0.04 


-0.00 


0.05 


1.2 - 1.4 


0.01 


0.03 


-0.00 


0.06 


1.4 - 1.6 


0.01 


0.03 


-0.01 


0.07 


1.0 - l.O 


a m 
U.Ul 


A A/) 

U.U4 


A A1 

-U.Ul 


A AA 

u.uy 


1.8 - 2.0 


0.01 


0.03 


-0.02 


0.11 


0.0 - 0.2 


0.01 


0.05 


0.01 


0.04 


0.2 - 0.4 


0.01 


0.05 


0.00 


0.04 


0.4 - 0.6 


0.02 


0.05 


0.00 


0.04 


0.6 - 0.8 


0.01 


0.05 


0.00 


0.05 


0.8 - 1.0 


0.01 


0.06 


0.00 


0.06 


1.0 - 1.2 


0.01 


0.05 


-0.01 


0.07 


1.2 - 1.4 


0.01 


0.05 


-0.01 


0.09 


1.4 - 1.6 


0.01 


0.05 


-0.01 


0.09 


l.O - l.O 


a m 
U.Ul 


A A^ 
U.UO 


A AQ 
-U.Uo 


U.lo 


1.8 - 2.0 


0.02 


0.05 


-0.03 


0.15 


0.0 - 0.2 


0.02 


0.08 


0.00 


0.05 


0.2 - 0.4 


0.01 


0.08 


0.01 


0.05 


0.4 - 0.6 


0.02 


0.09 


0.01 


0.06 


0.6 - 0.8 


0.02 


0.08 


0.00 


0.06 


0.8 - 1.0 


0.02 


0.06 


-0.00 


0.07 


1.0 - 1.2 


0.02 


0.07 


-0.01 


0.09 


1.2 - 1.4 


0.01 


0.08 


-0.01 


0.11 


1.4 - 1.6 


0.01 


0.08 


-0.02 


0.14 


1.0 - l.O 


a no 
U.Uz 


A AQ 

U.Uo 


A A/1 

-U.U4 


A 1 Q 

U. lo 


1.8 - 2.0 


0.01 


0.08 


-0.02 


0.18 


0.0 - 0.2 


0.02 


0.10 


0.02 


0.07 


0.2 - 0.4 


0.03 


0.12 


0.01 


0.08 


0.4 - 0.6 


0.03 


0.12 


0.01 


0.08 


0.6 - 0.8 


0.03 


0.11 


0.00 


0.09 


0.8 - 1.0 


0.02 


0.12 


0.01 


0.11 


1.0 - 1.2 


0.02 


0.11 


-0.02 


0.14 


1.2 - 1.4 


0.02 


0.13 


-0.01 


0.16 


1.4 - 1.6 


0.02 


0.13 


-0.02 


0.20 


i.o - i.o 


a aq 
U.Uo 


A 1 1 
U.ll 


A A/1 
-U.U4 


A Ol 
U. Zl 


1.8 - 2.0 


0.02 


0.12 


-0.07 


0.27 


0.0 - 0.2 


0.04 


0.15 


0.03 


0.10 


0.2 - 0.4 


0.03 


0.17 


0.03 


0.10 


0.4 - 0.6 


0.02 


0.18 


0.02 


0.11 


0.6 - 0.8 


0.02 


0.16 


0.01 


0.13 


0.8 - 1.0 


0.03 


0.14 


0.00 


0.14 


1.0 - 1.2 


0.03 


0.16 


-0.01 


0.19 


1.2 - 1.4 


0.01 


0.18 


-0.05 


0.24 


1.4 - 1.6 


0.03 


0.16 


-0.05 


0.26 


1.0 - l.O 


a no 
U.Uz 


U. 1 1 


A AC 

-U.Uo 


A QO 


1.8 - 2.0 


0.03 


0.16 


-0.10 


0.38 


0.0 - 0.2 


0.04 


0.22 


0.05 


0.13 


0.2 - 0.4 


0.06 


0.21 


0.03 


0.13 


0.4 - 0.6 


0.04 


0.22 


0.04 


0.15 


0.6 - 0.8 


0.05 


0.23 


0.01 


0.19 


0.8 - 1.0 


0.01 


0.25 


-0.01 


0.22 


1.0 - 1.2 


0.02 


0.24 


-0.03 


0.24 


1.2 - 1.4 


0.05 


0.24 


-0.05 


0.28 


1.4 - 1.6 


0.03 


0.24 


-0.15 


0.44 


1.6 - 1.8 


0.02 


0.20 


-0.11 


0.45 


1.8 - 2.0 


0.05 


0.21 


-0.13 


0.43 


0.0 - 0.2 


0.07 


0.26 


0.04 


0.15 


0.2 - 0.4 


0.10 


0.28 


0.05 


0.17 


0.4 - 0.6 


0.07 


0.26 


0.04 


0.20 


0.6 - 0.8 


0.08 


0.28 


0.03 


0.22 


0.8 - 1.0 


0.05 


0.28 


0.00 


0.27 


1.0 - 1.2 


0.07 


0.33 


-0.02 


0.33 


1.2 - 1.4 


0.08 


0.29 


-0.07 


0.37 


1.4 - 1.6 


0.04 


0.26 


-0.14 


0.48 


1.6 - 1.8 


0.06 


0.27 


-0.12 


0.48 


1.8 - 2.0 


0.06 


0.28 


-0.12 


0.49 



23.0 - 23.5.. 



23.5 - 24.0.. 



24.0 - 24.5.. 



24.5 - 25.0.. 



25.0 - 25.5.. 



25.5 - 26.0.. 



26.0 - 26.5.. 



a Recovered magnitude or color range of simulated galaxies. 

b Average value of the difference between measured and input magnitudes or colors. Significant departures from zero imply systematic biases in 
the photometry. 

c The standard deviation in the difference between measured and input magnitudes or colors. 
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Fig. 2. — Color-magnitude diagram for sources detected in the NB4980 image. Right axis indicates the rest-frame Lya equivalent 
width (assuming z = 3.09) implied by the BV-NB4980 color for sources at 3.054 < z < 3.12. Isophotes for W-NB4980 colors 
are defined in the NB4980 image, while NB4980 magnitudes are "Kron-like" (i.e., SExtractor MAG_AUTO) total magnitudes. 
The solid curve indicates our formal lower-limit threshold in BV— NB4980, though sources with measured values above this 
limit are plotted with their measured values, excepting those with measured colors BV — NB4980 > 3.75 which are plotted as 
BV — NB4980 = 3.75. Dotted lines indicate ±3a scatter in the data and dashed lines indicate our limits defining the LAE 
candidate sample. Circled points indicate sources with NB3640 detections. 
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26 LBGs in our sample qualify as LAEs by our criteria, 
though none of these are detected in NB3640. Informa- 
tion for LBGs and LAEs with no NB3640 detection is 
given in the appendix. 

Figure displays the smoothed NB3640, R and 
iKT/ACS-F814W images for the 6 LBGs with NB3640 
detections. Figures 0] - |S] display the smoothed NB3640, 
LyA, R and ffST/ACS-F814W images for the 28 (main 
and faint sample) LAEs with NB3640 detections. The 
contours in these Figures represent NB3640 flux levels 
corresponding to 28.81, 28.06 and 27.62 mag arcsec -2 (1- 
, 2- and 3-cr pixel-to-pixel fluctuations in the smoothed 
image). 

We also investigated the average i?-band and NB3640 
flux of various subsamples of LBGs and LAEs by stack- 
ing postage stamp images of the individual galaxies. 
Regions containing flux from sources obvious neighbors 
were removed using masks created from the SExtractor 
"segmentation" images. Both the i?-band and NB3640 
postage stamps were centered on the R detection for 
LBGs and on the NB4980 detection for LAEs and, after 
the rejection of masked pixels, averaged to create each 
stack. No significant NB3640 flux was detected in ei- 
ther the stack of the 20 LBGs (down to a 2tr limiting 
magnitude of 28.14) or of the 83 LAEs (down to a 2a 
limiting magnitude of 28.91) lacking individual NB3640 
detections^ The results from the photometry on stacked 
subsamples are reported in Table |H1 

As can be seen in Figures [3]- [8] (as well as Tables [4] and 
[5]), the NB3640, Lya, and _R-band fluxes are not precisely 
co-spatial in many of our sources. Spatial offsets between 
detected ionizing and non-ionizing UV continuum fluxes 
i.e., significant values of Ar), have been previous noted 
Iwata et alj|2009t llnoue et al.ll201lD . The values of A R 
in our LAE sample are relatively small (< 0'.'6), how- 
ever, compared to the values of Al v a and LyA-R offsets 
(Al v ar)- The median values for main sample LAEs 
with NB3640 detections are A R = 0'.'32, A LyA = 0'.'82 
and Al v a,r — 0'.'58. We show the individual values of 
Al v a versus A R , as well as A R versus NB3640 magni- 
tude, in Figure M 

3.5. The source of the NB3640 detections 

The detections in the NB3640 image coincident with 
the locations of z > 3.06 LBGs and LAEs are likely 
due to some combination of: (i) photons of wavelengths 
shortward of the Lyman limit escaping from the galax- 
ies and being redshifted above the Lyman break before 
the neutral phase of the IGM is able to absorb them (ii) 
non-ionizing flux from lower-redshift galaxies randomly 
coincident on the sky with the high-z systems. In princi- 
ple, the high spatial resolution of the HST/ACS imaging 
available for many of our sources with NB3640 detec- 
tions could help distinguish between the scenarios. How- 
ever, galaxies observe d in the rest-frame UV wavelengths 
often appear clum py (jLotz et al.l 120061 : ILaw et all 120071 : 
iPeter et al.|[2007h - i.e., as distinct regions significantly 
brighter than the regions surrounding them - particu- 
larly at high redshift. Thus, even with the resolution 

5 Although we did not test the NB3640 filter in the lab for the 
possibility of a "red leak" , the strong detection limits in the stacks 
indicate that there is no such leak at a level that could affect our 
results. 



afforded by HST/ ACS imaging, it is not generally possi- 
ble to distinguish in individual systems between clumpy 
emission from a single redshift (z ~ 3) and multiple 
sources at different redshifts. 

It is possible to address the issue of foreground contam- 
ination in a statistical manner, however. The positions 
of any foreground objects will be uncorrelated with that 
of the z > 3.06 sources, and their surface density can 
be predicted using, e.g., [/-band n umber counts in th e 
ultradeep VLT/VIMOS catalogue (|Nonino et all 120091 ). 
We predict the surface density of such sources, ps, from 
the iNonino et all (|2009l ) number counts by inter polating 
the values in Table 1 of iVanzella et all (|2010bl ) for our 
magnitude range 24.5 < U < 27.25. To define the posi- 
tions of the 26 LBGs, we use the centroids of the i?-band 
detections, which have confirmed spectroscopic redshifts 
z > 3.06. For the 110 LAEs, we first use the centroids 
of the detections in the LyA image (although, see be- 
low), which are very likely to be in the redshift range 
3.054 < z < 3.120. In Figure [10] we show the radial 
surface-density of all NB3640 detections surrounding the 
z > 3.06 sources, in bins of radial offset (i.e., for 
the LBGs and Al v a for the LAEs). The surface den- 
sity is defined as the number of detections divided by 
n{ r max- r min)N where r mm and r max are the minimum 
and maximum radii of the bin and N is the number of 
z > 3.06 sources, 26 for the LBGs and 110 for the LAEs. 
The left-hand panel shows the distribution for the LBG 
sample, while the bottom right panel shows the distribu- 
tion for the LAEs. The shaded portion of the histograms 
represent detections that we have identified as clearly 
belonging to a neighboring source other than the target 
LBG/LAE (and therefore do not appear in Tables [4] and 
[5]). The dashed and dotted lines indicate the predicted 
surface density of interlopers, ps, and the la uncertainty 
in the prediction, respectively. At large offset, the dis- 
tributions are comprised entirely of obvious neighbors, 
flat and consistent with ps- This consistency indicates 
that any possible difference between the J7-band num- 
ber counts and the surface density of NB3640 sources 
in our field (e.g., from the differences in filters, cosmic 
variance, etc.) can be ignored for our purposes. While 
disentangling true detections from interlopers becomes 
more difficult at smaller offsets, there is clearly a signif- 
icant excess of detections above the average foreground 
level, indicating that many of the retained NB3640 de- 
tections with offsets < 1" must be physically related to 
the redshift z > 3.06 sources. The magnitude of this ex- 
cess over that predicted by the [/-band number counts 
represents the predicted number of uncontaminated LyC 
detections, which we quantify below. 

We use ps to compute the full (radially-dependent) 
probability distribution for the number of int erlopers in 
our samples. Following Va nzella et al.l (|2010bl ). the prob- 
ability that a given LBG/LAE has a random foreground 
source with an offset in the range r m i n : r max is then 

Pr = K{r 2 max - r 2 mm ) X p s . (1) 

The probability that K of the N = 131 sources have a 
contaminant centered within the annulus is 

MK)=(%)p?(l (2) 
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TABLE 4 

Photometry for LBGs with NB3640 detections. 



ID 


RA 
(J2000) 


Dec 
(J2000) 




z h b 


H 


NB3640 


Ah c 


F UV d 
F LyC obs 


F UV e 

Ft „r 

cor 


MD46 


22:17:27.28 


0:18:09.7 


3.091 


3.080 


23.49 


25.22 


l'.'O 


4.9 ± 0.7 


2.4± 1.0 


C16 


22:17:31.95 


0:13:16.3 




3.0651 


23.62 


26.43 


1'.'9 


13.3 ±4.0 


6.8 ±3.3 


C49 


22:17:19.81 


0:18:18.8 


3.1629 


3.1492 


23.81 


26.84 


0'.'4 


16.4 ±6.1 


5.9 ±4.1 


D17 


22:17:18.86 


0:18:17.0 


3.0898 


3.0697 


24.29 


27.00 


0'.'9 


12.2 ± 5.1 


5.9 ±3.5 


aug96M16 


22:17:30.86 


0:13:10.8 


3.298 


3.285 


24.47 


25.23 


0'.'7 


2.0 ±0.4 


0.2 ±0.4 


MD32 


22:17:23.70 


0:16:01.6 


3.102 




25.14 


25.51 


0'.'4 


1.4 ±0.4 


0.6 ±0.3 



a Lya emission redshift. 

Instellar absorption redshift. 
c Spatial offset between the ccntroids of i?-band and NB3640 emission. 

Observed ratio and uncertainty in non-ionizing UV and LyC emission, inferred from the NB3640— R color. This value has not been corrected for 
either contamination by foreground sources or IGM absorption. 

c Ratio of non-ionizing UV continuum and LyC emission, corrected for both foreground contamination and IGM absorption. 



TABLE 5 

Photometry for LAEs with NB3640 detections. 



ID RA Dec 4980 W-NB4980 EW a R NB3640 A R b A LyA b p^- 





(J2000) 


(J2000) 






(A) 










'IsO obs 


LyC cor 


003 


22:17:24.79 


0:17:17.4 


23.20 


1.10 


39 


24.42 


24.74 


0'.'3 


o'.'e 


1.3 ±0.1 


0.6 ±0.3 


010 


22:17:20.38 


0:18:04.2 


24.03 


0.91 


28 


25.77 


26.74 


0'.'3 


0'.'3 


2.4 ± 1.1 


1.2 ±0.7 


016 


22:17:35.61 


0:18:00.2 


24.35 


1.66 


86 


26.24 


26.91 


0'.'9 


O'.'S 


1.8 ±0.9 


0.9 ±0.6 


018 


22:17:39.01 


0:17:26.4 


24.45 


1.53 


72 


26.25 


25.69 


O'.'l 


l'.'O 


0.6 ± 0.2 


0.3 ±0.2 


019 


22:17:26.15 


0:13:20.1 


24.47 


1.50 


70 


25.70 


26.24 


0'.'4 


0'.'9 


1.6 ±0.6 


0.8 ±0.4 


021 


22:17:18.77 


0:15:18.1 


24.57 


1.70 


92 


>27 


27.17 


1'.'3 


T.'4 


< 1.8 


< 0.9 


025 


22:17:36.74 


0:16:28.8 


24.69 


0.88 


27 


25.54 


25.85 


0'.'3 


T.'2 


1.3 ± 0.4 


0.6 ±0.3 


028 


22:17:31.80 


0:17:17.9 


24.77 


1.32 


54 


25.50 


26.71 


0'.'3 


0'.'9 


3.1 ± 1.3 


1.5 ±0.9 


034 


22:17:23.41 


0:16:35.4 


24.86 


0.89 


27 


25.42 


25.76 


0'.'5 


O'.'O 


1.4 ±0.4 


0.7 ±0.3 


038 


22:17:34.77 


0:15:41.3 


24.90 


0.84 


25 


26.17 


25.82 


O'.'l 


0'.'7 


0.7 ±0.3 


0.3 ±0.2 


039 


22:17:24.08 


0:11:31.7 


24.92 


1.03 


35 


26.48 


26.77 


0'.'4 


0'.'8 


1.3 ± 0.7 


0.6 ±0.4 


041 


22:17:24.54 


0:15:06.7 


24.93 


0.92 


29 


25.97 


25.94 


0'.'4 


0'.'5 


1.0 ±0.4 


0.5 ±0.3 


046 


22:17:21.47 


0:14:54.6 


25.04 


> 2.54 


>277 


>27 


26.43 


1'.'8 


1'.'8 


< 0.8 


< 0.4 


048 


22:17:27.37 


0:16:51.5 


25.06 


1.72 


94 


26.73 


26.00 


2'.' 2 


1'.'7 


0.5 ± 0.2 


0.2 ±0.2 


051 


22:17:33.72 


0:15:04.9 


25.10 


1.15 


41 


26.26 


27.21 


0'.'3 


0'.'7 


2.4 ± 1.3 


1.2 ±0.8 


053 


22:17:34.70 


0:16:33.4 


25.15 


0.96 


31 


26.53 


26.98 


0'.'6 


0'.'9 


1.5 ± 0.8 


0.7 ±0.5 


064 


22:17:35.42 


0:12:14.6 


25.38 


1.04 


35 


26.05 


26.74 


O'.'l 


l'.'O 


1.9 ± 0.9 


0.9 ±0.6 


069 


22:17:18.96 


0:11:12.0 


25.44 


0.71 


19 


24.62 


27.22 


0'.'2 


0'.'9 


10.9 ±5.1 


5.3 ±3.3 


074 


22:17:36.47 


0:12:54.8 


25.50 


0.75 


21 


26.16 


25.52 


O'.'l 


0'.'7 


0.6 ±0.2 


0.3 ±0.1 


077 


22:17:37.95 


0:11:01.3 


25.54 


0.73 


20 


26.01 


26.36 


0'.'2 


l'.'l 


1.4 ±0.6 


0.7 ±0.4 


081 


22:17:29.22 


0:14:48.7 


25.60 


1.45 


65 


>27 


26.79 


0'.'6° 


0'.'7 


< 1.2 


< 0.6 


083 


22:17:28.46 


0:12:08.9 


25.63 


1.48 


68 


26.46 


26.84 


0'.'4 


0'.'6 


1.4 ±0.7 


0.7 ±0.5 


084 


22:17:19.90 


0:15:14.9 


25.64 


0.71 


20 


>27 


26.50 


O'.'l 


0'.'8 


< 0.9 


< 0.4 


087 


22:17:37.07 


0:13:21.5 


25.72 


1.03 


35 


>27 


27.26 


O'.'l 6 


T.'6 


< 2.0 


< 1.0 


096 


22:17:38.93 


0:11:37.4 


25.83 


0.91 


28 


26.64 


26.45 


0'.'6 


0'.'6 


0.8 ±0.4 


0.4 ±0.3 


101 


22:17:25.33 


0:17:22.5 


25.90 


0.81 


24 


26.89 


26.58 


0'.'2 


T.'4 


0.8 ±0.4 


0.4 ±0.2 


102 


22:17:24.00 


0:16:27.6 


25.92 


0.78 


22 


25.91 


26.05 


O'.'l 


0'.'2 


1.1 ± 0.4 


0.6 ±0.3 


118 


22:17:35.28 


0:10:59.9 


26.18 


> 1.39 


>243 


26.20 


26.11 


0'.'4 


T.'4 


0.9 ±0.4 


0.4 ±0.3 



a Lya rest equivalent width estimated from BV— NB4980 color. 

b Spatial offset between the centroids of R or Lya and NB3640 emission. 

c Observed ratio and uncertainty in non-ionizing UV and LyC emission, inferred from the NB3640 — R color. This value has not been corrected for 
cither contamination by foreground sources or IGM absorption. 

Ratio of non-ionizing UV continuum and LyC emission, corrected for both foreground contamination and IGM absorption. 
c Offset determined from ccntroid of BV detection. 



The expectation value for the number of contaminated 
LBGs/LAEs is 

(n r )=Y,Kfr(K). (3) 

K 

Note that, dividing (n r ) by the areas corresponding to 
each bin in Figure [TU1 one recovers the surface density 
indicated by the dashed line. 
Having determined f r (K) for bins of offset, we can use 



Monte Carlo simulations to make a robust prediction of 
the number of uncontaminatcd NB3640 detections. In 
each bin, we randomly select a number of interlopers, 
Winter-, from f r (K). We then subtract the number of 
previously rejected sources in the bin, n re j (i.e., those 
corresponding to obvious neighbors, and indicated by the 
hatched histogram in Figure H0| . If N = rimter — n re j > 
0, we randomly reject N of the retained detections in 
the given bin. We repeated the simulation 1000 times 



10 



NESTOR ET AL. 



NB3640 



F814W 



aug96 
M16 



C16 



C49 



MD32 



MD46 



D17 



- / 



,1V, ^ 




1 J 



I 



i 



■■v, jj i " 

■ m 

* ' • ••• :• 



' ' M. 



Fig. 3.— LBGs with z > 3.06 detected in NB3640. Images are centered on the 7?-band centroid and span 7"x7". The NB3640 
filter (image shown here after smoothing) samples the rest-frame LyC, while the R and F814W filters sample the rest-frame 
non-ionizing UV continuum. 

to determine the average and la deviation of uncontam- 
inated NB3640 detections. The results of our simulation 
indicate that 3.2 ± 1.2 of the 26 LBGs and 18.0 ± 3.2 
of the 110 main sample LAEs have NB3640 detections 
that are uncontaminated by foreground sources. The two 
samples are not independent, as they share 5 sources; 
considering the sample as a whole, our simulation deter- 
mines that 21.6 ± 3.8 of the 131 sources have uncontam- 
inated NB3640 detections. These same simulations can 
also be used to contamination-correct the ensemble av- 
erage NB3640 magnitudes determined from the stacked 
images (Table [6]) or by summing the individual magni- 
tudes assuming zero flux from the non-detections (see 



sap- 

As discussed in §3.4[ detections in the NB3640 images 
associated with LAEs are less spatially consistent with 
the detected Lya emission than with the rest-frame UV 
continuum emission. Thus, it will be more suitable to de- 
fine the LAE positions using the centroids of the i?-band 
detections (i.e., use in place of Al v a) for the pur- 
pose of estimating the contamination, unless significant 
numbers of the NB3640 detections that we have associ- 
ated with z ~ 3.09 LAEs are interlopers and the cor- 
responding i?-band detection is dominated by the inter- 
loper. However, as we concluded above that the majority 
of such detections are not contaminated, the A# offsets 
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Fig. 4. — As in Figure but for LAEs with NB3640 detections. Also shown is the LyA (i.e. 
NB4980, see text) image. Images are centered on the LyA detection and span 7"x7". 



_Bl/-continuum subtracted 



should lead to more realistic predictions of the number of 
uncontaminated sources. The upper-right panel of Fig- 
ure[Tn]shows the radial surface density of detections using 
An offsets for the 78 LAEs for which an i?-band centroid 
was measurable, Abv offsets for 9 LAEs detected in the 
BV image but not R, and Anb4980 offsets for the re- 
maining 23 LAEs. Using these offsets, our Monte Carlo 
simulations predicts 20.8 ± 2.7 of the 110 main sample 
LAEs and 24.5 ± 3.0 of the 131 combined (LBG+LAE) 
sample sources are uncontamined. 

Finally, we consider the possibility that the LAEs 
with NB3640 detections having small NB3640-i? off- 
sets are not at z ~ 3.09, but rather, e.g., [Oil] emit- 
ters at z ~ 0.34. In this scenario, the NB3640 filter 



would sample the rest-frame UV flux and NB4980— BV 
would be a measure [Oil] emission, both of which should 
closely trace regions of active star-formation, while R 
would correspond to the rest-frame optical. The ob- 
served NB3640 flux would then be more spatially cor- 
related with NB4980— BV compared to R. As can be 
seen in Figure |H1 however, the opposite is the case. 

4. RESULTS 

One of the primary goals of this work is to estimate the 
comoving emissivity of ionizing radiation that escapes 
from star- forming galaxies at z ~ 3. Ideally, this esti- 
mate would entail deep, rest-frame LyC imaging of sev- 
eral large, independent fields to determine the A ~ 900A 
luminosity function; however, such an experiment would 
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Fig. 5. — Same as Figure [4] 




be prohibitively expensive in terms of observing time. 
Here, we have taken advantage of the large number of 
galaxies at z > 3.06 in the SSA22a field, which contains a 
known protocluster at z = 3.09, to measure the emergent 
UV flux escaping from a large sample of LBGs and LAEs 
both redward and blueward of the Lyman limit. The re- 
sulting ratios of non-ionizing to ionizing flux densities al- 
low us to use the previously measured global non-ionizing 
UV emissivity at z ~ 3 to determine the overall contri- 
bution to the ionizing UV background at these redshifts 
from galaxies represented by those in our subsamples. 



4.1. Observed UV to Lyman- continuum flux density 

ratios 



The Suprime-Cam R filter has an effective wavelength 
of - 1600A and FWHM of - 300A in the rest-frame 
of our LBG and LAE samples. This is is a useful spec- 
tral range for determining UV to LyC flux density ratios, 
as the luminosity functi on at A ~ 1700 A h as been well 
studied at z ~ 3 (e.g., iReddv et all 12008ft . We show 
the NB3640— R colors of our sources as a function of 
R in Figure [TT] Also shown are the colors of stacks 
of all LBGs, all LAEs, LBGs and LAEs having 
no individual NB3640 detection, and subsamples LAEs 
in bins of R. Sources with NB3640 detections in ob- 
jects brighter than R ~ 25, while few (5 LBGs and 2 
LAEs), span a large range NB3640-i? color (~ - 3). 
While the apparent trend in color with R is in part a 
result of the NB3640 detection limit (indicated by the 
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Fig. 6. — Same as Figure |U 




dotted line in Figure [TTj) and the increased photometric 
uncertainty for fainter systems (e.g., ONB3640-R ^ 0.5 for 
R ~ 26.5, NB3640-i? ~ 0; see Table©, photometry on 
the stacks of LBGs and LAEs in bins of R indicates a 
significant difference in the average NB3640— R color of 
systems brighter and fainter than R ~ 25. 

The observed UV to LyC flux density ratios can be 
computed directly from the values of NB3640 and R. 
We present these ratios in Tables 2] and [5j Our value for 
C49, Fuv/FhyC — 16.6 ±6.1, is c onsistent with that de - 
rived from the LRIS spectrum bv lShaplev et al. | ( 2006). 
Fuv / FlvC = 12.7 ± 1.8. However, unlike IShaplev et all 
(2006) we do not detect D3 in NB3640, consistent with 
the findings of llwata et all (|2009f ). The average ratios 
for the ensembles of LBGs and LAEs are also of inter- 



est. These can be computed either from the photome- 
try on the stacked images or from the individual R and 
NB3640 measurements. Motivated by the non-detection 
of NB3640 flux in the stacks of individually non-detected 
sources, we assume zero flux for the non-detections in 
NB3640 when computing the ensemble ratios from in- 
dividual measurements!.] Doing so, we find values of 
{Fuv I F LyC ) obs = 17.8i^ 6 4 and 4.2±i;g for the ensem- 
bles of LBGs and main sample LAEs, respectively. Con- 
sidering only sources with NB3640 detections, we find 

{FuvlF LyC ) d Z = 5.2±i;2 for LBGs and 1.3 ± 0.3 for 
LAEs. In all cases, the errors are dominated by the sam- 



6 For ease we also assume zero flux in R for LAEs with no _R-band 
detection. This assumption has negligible effect on our results. 
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Fig. 7. — Same as Figure [4] 




pie variance which is determined using bootstrap resam- 
pling. The ensemble results are summarized in Table [3 
The ratios determined from the stacked data, which we 
report in Table |H1 are consistent with the above values 
but with larger uncertainties. For the sources with non- 
detections, we find limits of {Fuv / FL y c) r ^~ det > 33.7 
and > 12.2 for the LBGs and LAEs respectively. 

4.2. IGM absorption of Lyman- continuum flux 

Any LyC flux that escapes from the galaxies in our 
samples is diminished by absorption from the neutral 
phase of the IGM. In order to model IGM opacity, we 
use Monte Carlo simulations of intergalactic absorption 
to generate large (500) samples of random sightlines to 
the redshifts of our target galaxies. Our simulations are 



a mod ified version of those presented in IShaplev et al.l 
((200^ . 

Absorbers are drawn at random from column-density 
and redshift distributions consis tent with recent de- 
terminations in the l iterat u re (llnoue fc Iwatal 120081 : 
iFaucher-Giguere et al.l 120081: ISongaila fc Cowid I2010D 
for both the lower-column-density Lya forest, and 
higher-column-density Lyman limit systems (LLSs) and 
Damped Lya systems (DLAs). Of particular rele- 
vance for estimating the attenuation just below 912 
A are the statistics for LLSs (N m > 10 172 ). For 
these systems, we assume a power-law distribution in 
column-density, proportiona l to -/V^j 1 ' 3 (|Inoue fc Iwatal 
20081: iProchaska et al.l [2010D . redshift evolution propor- 
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TABLE 6 
Photometry in stacked-images. 



sample 


num. 


R 


NB3640 a 


phot b 
cr NB3G40 — R 


sample c 
CT NB364()-fl 


NB3640-i? 


(F uv /F LyC ) obs d 


LBGs 


26 


24.23 


27.94 


0.41 


0.91 


3.71 ± 1.00 


30-5±S:i 


LBGs, non-detect 


20 


24.32 


> 28.14 






> 3.82 


> 33.7 


LBGs, detect 


6 


23.98 


25.81 


0.13 


0.23 


1.83 ±0.26 


■5 4 +i.5 


LAEs 


110 


26.14 


27.49 


0.14 


0.32 


1.35 ±0.35 


o cr+1.3 
°'°-1.0 


LAEs, non-detect 


83 


26.19 


> 28.91 






> 2.72 


> 12.2 


LAEs, detect 


27 


26.04 


26.22 


0.14 


0.57 


0.18 ±0.59 


1 2+ 1 ' 9 


LAEs, R < 25 


8 


24.44 


27.16 


0.36 


0.94 


2.72 ± 1.01 


12.2+7 8 4 s 


LAEs, 25 < R < 26 


23 


25.64 


26.69 


0.15 


0.26 


1.05 ±0.30 


2 6+ ' 8 
z - D -0.6 


LAEs, 26 < R < 27 


43 


26.68 


27.59 


0.27 


0.53 


0.91 ±0.59 


2 3+ 17 


LAEs, ij > 27 


36 


27.53 


28.41 


0.67 


0.84 


0.88 ± 1.07 


if (i+3.8 
Z - Z -1A 


faint LAEs 


20 


27.93 


> 28.13 






> 0.20 


> 1.2 



a Limits, when quoted, are 2a. 

Photometric color uncertainty determined from stacks of randomly sampled blank regions of the R and NB3640 images. 
c Sample color uncertainty determined by bootstrap resampling of the subsample. 

Observed ratio and uncertainty in non-ionizing UV and LyC emission, inferred from the NB3640— R color of subsample stacks. This value has 
not been corrected for either contamination by foreground sources or IGM absorption. 
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1 2 

NB3640-R offset (arcsec) 

Fig. 9. — Left: Comparison of the spatial offset between detections in the NB3640 and LyA images with those between NB3640 
and R for our LAE sample. The smaller NB3640-i? offsets relative to NB3640-Ly J 4 strongly favors the sources being LAE 
emitters at z ~ 3.09 over, e.g., [Off] emitters at lower redshift (see text). Right: NB3640-/? spatial offsets versus NB3640 
magnitude. Open points represent LBGs while solid points represent LAEs. 
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Fig. 10. — The radial surface- density of NB3640 detections around galaxies in our LBG and LAE samples. The solid histograms 
includes all sources identified by SExtractor; the subset of these sources associated with obvious neighbors (which we have 
excluded from our catalog, see t|3.3p are represented by the hatched region. The dashed lines indicates the global surface density 
of sources in our NB3640 magnitude range and thus represents the expected level of contamination, while the dotted lines 
represent the expected la scatter in the contamination. The excess surface density at low offsets indicates that several of our 
low offset LBG NB3640 detections and the bulk of the low offset LAE NB3640 detections are physically associated with the 
z ~ 3.09 sources and not random foreground interlopers. The top panels use the displacement of the LyC centroid from the that 
of the J?-band detection (i.e., Ar) for LBGs and LAEs with ii-band detections, the displacements from the BV centroid (Abv) 
for LAEs detected in BV but not R, and from the NB4980 (Anb498o) detection for the remaining LAEs. The bottom panel 
uses displacements from the LyA detections (Al v a) of LAEs. The excess surface density is contained within a larger range of 
offsets in this panel, as the An and Abv values tend to be smaller than the Al v a values in our LAE sample (see Figure |9}. 
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Fig. 11. — Rest-frame LyC— UV color versus rest-frame UV magnitude. Open (filled) points and arrows represent LBGs (LAEs 
with R and/or NB3640 detections). The uncertainty in color for fainter systems is fairly large; e.g., <7nb3640-.r > 0.5 for R ~ 26.5, 
NB3640— R ~ (see Table [2]). Large points represent measurements made on stacks of the entire sample of LBGs and LAEs. 
Lower-limits are for stacks of LBGs (far-left arrow) and LAEs (middle arrow) with no individual NB3640 detection. Smaller 
error-bars indicate measurement uncertainty, while the larger bars also include sample variance computed through bootstrap 
resampling. Red crosses (and far-right lower-limit) are measurements from stacks in bins of R magnitude. Note that five of the 
eight LAEs in the brightest bin are also LBGs. The top axis is the corresponding absolute magnitude scale for z = 3.09. 



TABLE 7 

Average UV to LyC flux density ratios. 



Correction 


(NB3640) - (R) a 


LBGs 

(Fuv/F Ly o) b 


LAEs 

(NB3640) - (R) a 


{Fuv/Flvc) b 






Full Ensembles 






none 

contamination c 

IGM + contamination d 


3.13 ±0.50 
3.68 ±0.69 
2.64 ±0.70 


17 O + 10.4 
„ q+10.3 


1.56 ±0.28 
1.82 ±0.36 
0.86 ±0.36 


4 2+ 1 ' 2 

r; 4+2.1 

9 9 +0.9 
z ' z -0.6 



Sources with NB3640 detections only 



none 

contamination c 

IGM + contamination d 



1.78 ±0.33 
1.65 ±0.51 
0.45 ± 0.59 



5.2 
4.6 
1.5 



1.8 
1.4 
+2.7 
1.7 
1.1 
0.6 



0.25 ±0.25 
0.28 ±0.30 
-0.69 ±0.32 



1.3 ±0.3 

i ,+0.4 
L - y -0.3 
5+ - 2 



a Color determined from average NB3640 and i?-band fluxes. Uncertainties include individual flux and sample uncertainties. 

Ratio and uncertainty in non-ionizing UV and LyC emission inferred from (NB3640) — (R) color. 
c Color and flux density ratio after statistically correcting sample for foreground contamination of NB3640 fluxes. 

Color and flux density ratio after correcting sample for forground contamination and IGM absorption of NB3640 fluxes. 
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tional to (1 + z) 2M (|Songaila fc Cowid 12010ft . with a 
tot al number per un it reds hift at z ~ 3 of dNj dz « 
2 (ISongaila fc Cowid [20101: iStengler-Larrea et all [19951: 
|Steide]||1992l). Randomly generated absorbers are then 
applied to the unabsorbed continuum of each model 
sightline until the total number of absorbers is equal to a 
random Poisson realization of the integral of the absorber 
column-density and redshift distributions over the rele- 
vant ranges in column-density (-/Vhi = 10 12 — 10 22 cm~ 2 ) 
and redshift (z = 1.7 — z source , where z source is the 
redshift of the target galaxy). The simulations pre- 
dict a mean free path to ionizing radiation at z ~ 3 
(as estimated from the redshift at which F cont /Fgi2 = 
exp(— 1), following Proch aska et al.l (|2009l )L which is 
consistent with recent deter minations in the literature of 
A m fr, ~ 72—85 proper Mpc (iFaucher-Giguere et al.ll2008t 
iProchaska et al.l 120091: ISongaila fe Cowid 12010ft . Using 
the relation given bv ISongaila fc Cow ic (2010), we adopt 
a value of X m fp = 75.6 Mpc at z = 3.1 for the analysis 
in EH 

For each target redshift, the resulting sample of 500 
model sightlines is analyzed to obtain the median and 
standard deviation of the attenuation factor within the 
NB3640 filter. This quantity differs from the one pre- 
sented in lShaplev et al.l (|2006l ). which corresponds to the 
attenuation factor in the fixed rest-frame bandpass of 
880 - 910 A. In iShaplev et al.l (|2006l ). a uniform trans- 
mission function was also assumed across this fixed rest- 
frame bandpass. Here we estimate the attenuation in a 
fixed observed frame bandpass, and also take into account 
the shape of the NB3640 filter. For higher redshifts, the 
attenuation within the NB3640 filter is stronger not only 
due to the the evolving physical properties of the Lya for- 
est, but also the fact that the NB3640 filter probes bluer 
rest-frame wavelengths sensitive to photons that have 
traveled a longer path-length from the emitting galaxy. 

For the six LBGs with NB3640 detections, the median 
and standard deviation of the 500 attenuation factors de- 
termined for the appropriate redshift are used to correct 
the observed UV to LyC flux density ratios. We assume 
z = 3.09 for each LAE. We use the average redshift of 
the LBG sample, (z) = 3.10, and of the subsample with 
NB3640 detections, (z) — 3.13, to correct the subsample 
ratios. In Figure IT21 we show the distribution of attenu- 
ation factors for our z — 3.09 simulation as well as the 
distributions for the LBG and LAE ensemble average at- 
tenuation factors. When computing the uncertainties in 
the median attenuations for LBG and LAE subsamples, 
we assume the sightlines to the individual sources are 
uncorrelated. The corrected flux density ratios are pre- 
sented in Tables [4] and [5] for individual sources, and in 
Table [7] for ensembles. 

4.3. The escape fraction of ionizing photons 

The UV to LyC flux density ratios determined in our 
samples of LBGs and LAEs can be used to estimate 
the fraction of LyC photons that escape from galaxies 
into the IGM, f^P '. However, we must first correct the 
ensemble average flux density ratios for contamination 
( §3.5j) and IGM absorption ( §4.2|) . In each realization 
of the Monte Carlo simulations described in M3 . 5 1 we set 
the NB3640 fluxes to zero for sources selected as inter- 
lopers by the simulation. We retain all measured i?-band 
fluxes, which is an appropriate approximation if the in- 
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Fig. 12. — Results of Monte Carlo simulations of IGM opac- 
ity along the line of sight to our high-redshift targets. To ac- 
count for IGM opacity, we generated a large number of ran- 
dom realizations of our measurements through the NB3640 
filter. These simulations were tuned to the observation of an 
individual sightline (relevant for observations of individual 
LBGs and LAEs) as well as the average of ensembles of sight- 
lines (relevant for sample averages and stacked images). In 
each simulation, using 500 random realizations, we recovered 
the distribution of the ratio between observed and unatten- 
uated NB3640 flux, denoted here as "NB3640 Attenuation 
Factor." In each panel, the normalized distribution of sim- 
ulated NB3640 attenuation factors is plotted, relevant for a 
given observed quantity. The dashed vertical line indicates 
the median NB3640 attenuation factor from each simulation, 
used to correct the observed individual or ensemble average 
Fuv I Fl V c flux-density ratio for IGM attenuation. (Top) The 
results for a single sightline to a source at z = 3.09. (Middle) 
The results for the ensemble average of 26 LBGs, with a mean 
redshift of (z) = 3.10. (Bottom) The results for the ensemble 
average of 110 LAEs, with a mean redshift of (z) — 3.09. 

terlopers are significantly bluer in observed NB640— R 
relative to the z ~ 3 source, and in any case is the most 
conservative approximation. We then apply the appro- 
priate IGM correction ( tJ4.2[) to the average flux density 
ratios determined from the contamination-corrected sam- 
ples. We present the raw and corrected colors and flux 
density ratios for the ensembles and subsamples with 
NB3640 detections in Table [71 Errors include sample 
variance computed by first randomly reassigning indi- 
vidual magnitudes based on the measured magnitude and 
error, assuming Gaussian magnitude uncertainties deter- 
mined from our photometric simulations f £|3.ip . and then 
bootstrap resampling the data set. After applying both 

corrections, we find (F uv / 'F Ly c) JL = 11.3±l°f for our 
sample of LBGs. Our sample of LAEs has on average 
much bluer LyC— UV colors, with a corrected flux den- 
sity ratio of {F uv / F LyC ) L J r r = 2.2±°J. The corrected 
ratios represent the average relative flux density in the 
vicinity of z ~ 3 galaxies. 

The absolute LyC escape fraction can be inferred from 
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the relative escape fraction as f^ c C = f r e f c x f^V , where 
is the escape fraction of non-ionizing UV photons 
and flf c is determined from the corrected flux density 
ratio: 

rrel _ ( L Uv/ L Lyc) intr ( 
J esc rjp I T-i \ > V 1 ) 

and thus 

feJf = (Fuv/F Ly c) corr x (L uv / L LyC ) intr x j^v ^ 

where (Luv / 'L Lyc) - nt is the (unknown) intrinsic ratio 
of UV to LyC l uminosity densities produced in star- 
forming regions. IReddv et al.l (|2008[) report an average 
UV escape fraction of ~ 20 - 25% in z ~ 3 LBGs. 
While direct measurements of (Luv I L Lyc) ■ t are lack- 
ing, spectral synthesis models predict values of ~ 6 for 
reasonable assumption of stellar population ages, m etal- 
licities and initial mass functions (|Siana et~^l200l . We 
adopt these numbers as our fiducial values, resulting in 
expressions for the relative and absolute escape fractions: 
flf c = 6.0 x (Fuv/FLyc)-^ and tfjf = 0.225 x f£. 
We caution, however, that sources with escaping LyC 
flux may also have larger UV escape fractions. In any 
case, both f^V and (L uv /L LyC ) intv are highly uncer- 
tain and likely vary from object to object. Also note that 
our fiducial value for (Ljjy / Li,yc ) irifr differs fro m that 
used by iSteidel et afl (|2001ft and iShaplev et all (|2006l ). 
who adopt a value of three. 

Most of the 20 LBGs with non-detections in NB3640 
have only weak individual upper limits on the relative 
escape fraction, ranging from f r e f c < 0.5 to ^unity. The 
(2<t) upper limit on their average escape fractions are 

{fl?cY° n ~ det < 0.5 and {f^T ^ < 0.10. Of 
the six LBGs with NB3640 detections, three have cor- 
rected UV to LyC flux density ratios consistent with 
fete — However, the other three and the subsam- 
ple average have NB3640— R colors suggesting values 
greater than unity and therefore in conflict with current 
models. If we assume that LBGs with NB3640 detec- 
tions preferentially sample clear sightlines through the 
IGM and we can thus ignore the IGM correction (but 
retain the statistical contamination correction) we find 

(fefe) det > 0-82 (la) and (f^ v c c ) det > 0.18 for the sub- 
sample of LBGs with NB3640 detections. For the LBG 
ensemble (detections and non-detections) average, the 

value of (Fuv/FLyc)^ = 11.3^° 4 3 obtained by apply- 
ing both IGM absorption and contamination corrections 
implies average escape fractions of (/,T s e c) = 0-53^0 25 an d 
(fLvp) = O.ntoll- Note that all of the above uncer- 
tainties do not include contributions from the uncertain- 
ties in (L uv /L Lv c) mtr or . 

Due to the overall faintness of the LAE sample, we 
can place a meaningful lower limit on flf c for only a few 
of the individual systems having no NB3640 detection. 
For the non-detection subsample average, we determine a 
weak upper limit of (f^tc) ^ 1-2 for our fiducial intrinsic 
Lyman break of six. The LAEs with detections have ex- 
tremely blue NB3640— R colors, however, implying indi- 
vidual flf c values greater than unity at high significance 
for all but one system, even if we neglect corrections 
for IGM absorption. These extremely blue colors lead 
to an inferred ensemble relative escape fraction that is 



greater than unity at ~ 4<r. However, an intrinsic value of 
(L uv /L LyC ) intr = 6 may be high for the LAEs at z ~ 3, 
as their UV lum inosity appears to b e dominated by par- 
ticularly young ([Gawiser et al.ll2007l) and low-metallicity 
(jOno et al.M2010D populations. Thus, we also consider a 
lower intrinsic flux density ratio of three, which leads to 
{fife) > 0- 75 ( 2cr )> and {fel C ) > 0-17 for the ensemble. 
We note, however, th at as LAEs have bee n found to ex- 
hibit little reddening (jGawiser et al.ll2007t) . is likely 
larger in our LAE sample than the fiducial LBG value 
of 20 - 25%; thus our determinations of f^ s v c c for LAEs 
should be considered conservative lower-limits. Nonethe- 
less, the NB3640-.R colors of LAEs with NB3640 detec- 
tions are uncomfortably blue, even when neglecting the 
IGM-absorption corrections. The reas ons for such small 
observed flu x density ratios ( see also llwata et"aLl 120091 : 
linouil20T0l : llnoue et al.l 120111) in our LAE samples are 
not clear. 

4.4. The space density of ionizing photons at z ~ 3.09 

Based on the average corrected UV to LyC flux density 
ratios for ensembles of LBGs and LAEs f ij4.3[) . we can 
compute the comoving density of ionizing photons that 
escape from galaxies z ~ 3, e LyC = (F uv / 'F Ly c)~ orr euv, 
where euv is the comoving luminosity density in the non- 
ionizing UV continuum. The value of euv can be ob- 
tained by integrating the rest-frame UV luminosity func- 
tion. Thus, 

e Ly c = (Fuv/F LyC )^ rr 4>* L* T(2 + a, L mm /L*) (6) 

where T is the incomplete gamma function, L* is the lu- 
minosity corresponding to the characteristic magnitude 
M*, and L m i n is the minimum luminosity to which the 
LF is inte grated. The A ~ 1 700 A LF has been deter- 
mined by IReddv et all (|2008f) for LBGs at z ~ 3 with 
Schechter function parameters <fi* = 1.66 x 10~ 3 Mpc -3 , 
M* AB = -20.84, a = -1.57 and an effective redshift 

Z e ff = 3.05. 

The average emergent UV to LyC flux density ratio 
for our sample of LBGs was determined from galax- 
ies brighter than the photometric limit for spectroscopic 
confirmation of LBGs in the SSA22a field, R = 25.5. 
This limit corresponds to ~ 0.5L* at z ~ 3. Thus we 
adopt (F uv I F L V c)cfrr = ll -3-l°4 for galaxies with A ~ 
1700A continuum magnitudes L > 0.5L*. The galaxies 
in our LAE sample are generally fainter, reaching at least 
i?~27(~0.1i*). However, LAEs comprise only a sub- 
set of faint z ~ 3 galaxies and may provide a biased 
view of the average far-UV properties of galaxies with 
A ~ 1700A continuum magnitudes of 0.1L* < L < 0.5L*. 
We therefore estimate the contribution to thyC from 
fainter (O.li* < L < 0.5L*) galaxies using in turn both 
the flux density ratio determined from our LBG sam- 

pie ((F uv /F LyC )t B rr = n - 3 -5 4 3 ) and that determined 
from our LAE sample ((F uv /F LyC )^ = 2.2±°;°). Us- 
ing the LBG-determined flux density ratio, we compute 
eLyc(L > 0.5L*) = 7.7±3 7 x 10 24 erg s" 1 Hz" 1 Mpc" 3 
and e LyC (0.lL* < L < 0.5L*) = 9.2+^ 3 2 8 x 
10 24 erg s _1 Hz" 1 Mpc -3 , for a total of e LyC (L > 
OIL*) = le^g^xlO 24 ergs^ 1 Hz^ 1 Mpc" 3 . Consider- 
ing instead the ratio determined from our LAE sample for 
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the fainter galaxies, we find CL y c(0.1L* < L < 0.5L*) = 
47.0±il;f x 10 24 erg s" 1 Hz" 1 Mpc" 3 , for a total of 

£Lyc(L > 0.1L*) = 54.7tu°2 x 10 24 erg s" 1 Hz" 1 Mpc" 3 . 

At z ~ 3, estimates of the UV to LyC flux density ratio 
in star -forming galaxies have been made bv lSteidel et al.l 
(120011) using a comp osite spectrum of 29 LBGs and by 
Shapl ev et al.l (|2006f) from spectra of 14 LBGs in the 
SSA22a field. We show our estimates of tLyc(L > 
0.1L*) in Figure [131 together with estimate s that we de- 
rive u sing the flux density ratio s given by iSteidel et al.l 
([20f)l and lShaplev eiTaTI IpOOfih and assuming the same 
non-ionizing UV luminosity density from L > 0.1L* 
galaxi es that we computed above using the IReddv et al.1 
(2008) LF. Additionall y, using HST far-UV ima ging of 
galaxies at z ~ 1.3, ISiana et al.l (|2007l . I2010D deter- 
mined (3d) lower limits of (Fuv / Fl v c) corr > 168 and 
(F uv I F LyC ) corr > 63 on the average non- ionizing UV to 
~700A flux density ratio for 15 L ~ L* starburst galaxies 
in the GOODS fields and 21 primarily L < L* galaxies 
in the HDF-N and HUDF, respectively. We also show 
the upper-limit on eu v c(L > 0.1L*) that we derive from 
their flux ratio limits in Figure [T3l using the A ~ 1500A 
LF pa rameters determined at z ~ 1 from lArnouts et al.l 
( 2005 ) and integrating down to L = 0.1L*. The con- 
tribution to e r. v n from broad-line QSOs determined by 
iCowie et~a l. (2009) is also shown, for comparison. 

We can use our estimates of thyC to calculate the cor- 
responding (proper) metagalactic specific intensity, J„: 



Jv 



(1 + zft LyC 



47T 



(7) 



0.33±°;» x 



We find J V (L > 0.5L*) 
10~ 21 erg s _1 Hz -1 cm -2 sr" 1 and, assuming the LBG 
flux density ratio holds down to L — 0.1L*, J V (L > 
0.1L*) = 0-71±g;|| x 10~ 21 erg s" 1 Hz" 1 cm" 2 sr" 1 . 
Alternatively, computing the contribution from fainter 
galaxies assuming the flux density ratio determined 
from our LAE sample leads to J U (L > 0.1L*) = 
2.30±g |g x 10~ 21 erg s" 1 Hz" 1 cm" 2 sr" 1 . 

Furthermore, if we assume a form for the spectrum of 
escaping ionizing flux, we can estimate the contribution 
to the hydrogen photoionization rate Thi in the IGM 
from star- forming galaxies at z ~ 3. The photoioniza- 
tion rate is of particular interest as Thi has been in- 
ferred in an independent manner from Lycn forest dat a 
(e.g., IBolton et all [2001 iFaucher-Giguere et al.l l2008h . 
If we assume a power-law spectrum J„ = J„ (i//Vo) Qs , 
where vq corresponds to the Lyman limit and cthi = 
6.3 x 10 -18 cm -2 is the atomic hydrogen photoioniza- 
tion cross section at vq, then 



HI 



4 7T cr H I Ju 

h(3 — a s ) 



(8) 



where h is Planck's constant. The intrinsic spec- 
trum of star-forming galaxies below the Lyman limit is 
poorly constrained. Typical v alues adopted in past work 
have ranged from a s — —3 (IBolton fe Haehneltl l2007t 
Ouchi et all 120091 ) to a s = +0.5 ([Faucher-Giguere et al l 
2008). If we assume a value of a s — —3 we find a con- 
tribution of T H i = 0.6±o j x 10~ 12 s- 1 from L > 0.5L* 
galaxies. For the full range of L > 0.1L*, we deter- 



mine r H i = l-4±J;f x 10" 12 s- 1 using the LBG flux 
density ratio, or, computing the contribution for faint 
galaxies using the LAE ratio, r H i = 4.5^^ x 10~ 12 s" 1 . 
Estimates of Thi from Lyq forest optical depths ra nge 
from ~ 0.6 x 10" 12 s " 1 (|Faucher-Giguere et al.ll2008[ ) to 
~ 1.3 x 10" 12 s" 1 (jKirkman et~all 120051 ) . The value 
obtained by applying the LBG flux density ratio over 
L > 0.1L* is consistent with this range, while attribut- 
ing the LAE flux density ratio to fainter galaxies re- 
sults in an estimate that is 3<r larger than implied by 
the Lya forest data. In Figure [13] we show the val- 
ue s of ej,„c corresponding to the determinations of Th t 
bvlBolton fc Haehneltl (|2007t l and lFaucher-Giguere et al.1 
(2008) from z=2 to z ~ 5. We arrive at these values us- 
ing equations [7] and [S] together with our adopted values 
of as and A m / P . Our estimates of eL V c, Jv, and Thi are 
summarized in in Table [8j 

Finally, it should be noted that the values estimated in 
this subsection are sensitive to the adopted parameters. 
Specifically, integrating the LF down to L m i n — would 
result in a factor of two increase in tLyCi Jv, an< l Thi 
compared to an integ ration lower-limit of L^ m = 0.LL*. 
Using the UV LF of IReddv fc Steidell (l2009h would not 
change our estimates of the contribution to tLyC , Jv, 
and Thi from galaxies with L > 0.5L*, but would in- 
crease our estimates for the ranges L > 0.1 and 0.1L* < 
L < 0.5L* by 14% and 27%, respectively. Contrastingly, 
using the LyA offsets to determine the contamination 
correction for the LAEs fi )3.5p would decrease our esti- 
mates of the LAE contribution to tLyC , Jv, and Thi 
by ~ 15%, while adopting a far-UV spectral slope of 
a s — +0.5 would increase our estimates of Thi by a fac- 
tor of 2.4. 

5. DISCUSSION 
5.1. Summary of the LBG/LAE subsamples 
5.1.1. LBGs 

Of the 26 LBGs at z > 3.06 in the field of our NB3640 
image, we have detected six down to a limiting mag- 
nitude of 27.3. Correcting for foreground contamina- 
tion and absorption by the neutral IGM, we determine 

(F uv / F Ly c)^ orr = H-3t5°4 3 . Previous direct measure- 
ments of the relative UV to LyC emergent flux den- 
sity r atios at z ~ 3 have found values ra nging from 
^ 4 (ISteidel et all f200lL llwata et al.l [2001 ) to ~ 20 
(jShaplev et all 120061 ). Our value is roughly consistent 
with both extremes. For an intrinsic ratio of six and a 
non-ionizing UV escape fraction of ~ 20 - 25%, this range 
in (Fuv I Fl v c) corr corresponds to a range of 0.07 - 0.4 
in (ftyc), with our value giving (f^ c c ) = 0.12?%&. 

Monte Carlo simulations based on the distribution of 
offsets between the NB3640 and i?-band detections indi- 
cate that two to four of the LBGs with NB3640 detec- 
tions are uncontaminated by foreground sources. Thus 
our corrected LBG detection rate is 8 - 15% (la). Three 
of the detections have ratios of non-ionizing to ionizing 
flux densities that are in tension with that allowed by 
stellar population models. If we take the most conser- 
vative interpretation of our results and consider only 
C49 and D17 - which have A R < 1" and flf c < 1 - 
as uncontaminated detections, we arrive at an ensemble 

(fefc) LBG - 0-10 and (f^ v c c ) LBG a 0.02. However, our 
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Fig. 13. — Contributions to the comoving emissivity of ionizing radiation at redshifts < z < 5. The open square is our 
determination of the contribution to th y c from galaxies with L > 0.1L* assuming the estimate of the LBG flux density ratio 
holds over this luminosity range, while the filled square adopts our estimate of the LAE ratio for galaxies with Q.1L* < L < 0.5L*. 
Also sh own are estimates bas ed on the UV to LyC flux density ratios determined f rom: (star) spectra o f 14 LBGs in the SSA22a 
field bv lShaplev eTall (HI); (t riangle) the composite spectrum of 29 LBGs by ISteidel et all (|2001l ); and (upper limit) rest- 
frame ~ 700A HST images of 36 s tarburst galaxies at z « 1.3 bv lSiana eFaLl (|2007l . [2010l ). The red crosses are the contribution 
from broad line QSOs taken from lCowie et ahl l|2009l ). The open and filled circles are estimates of the total ionizing emissivity 
determined using equations [7] and [5] a nd estimates of the hydrogen photoionization rate in the Lya forest, respectively, by 
iFaucher-Gieuere et all l|2008f ) and lBolton fc Haehneltl (|2007l ). Upper limits are 3er. 



TABLE 8 

Contributions to the ionizing backgrounds. 



(Fuv/FLyc) a 


Luminosity range b 


CLyC C 

(xl0 24 ergs s" 1 Hz" 1 Mpc -3 ) 


I d 

(xl0~ 21 ergs s _1 Hz -1 cm -2 sr _1 ) 


Thi c 

(xlO-^s- 1 ) 


11 n + 10.3 

AJ - a -5.4 
9 9 +0.9 
z - z -0.6 


L/L* > 0.5 
0.1 < L/L* < 0.5 


7 7+7.1 
'■'-3.7 


n nn+0.30 
U.OO_g 16 

1 O7+0-74 
k - y < -0.57 


fi+ ' 6 

3.9±i;5 




Total 


54.7+";°, 


n m+0.80 
Z - OU -0.59 


4 5+i.e 
l.i 


11 n + 10.3 

11 " 3 -5.4 


L/L* > 0.1 


1G.911 5 , 5 


u - (± -0.34 


1 4+1-3 



a Corrected flux density ratio assumed for the luminosity range. 

Luminosity range over which the LF is integrated. 
c Comoving specific emissivity of ionizing radiation. 
^ Proper metagalactic specific intensity of ionizing radiation. 

c Inferred contribution to the intcrgalactic hydrogen photoionization rate from galaxies in the luminosity range. 



contamination simulations suggest that this scenario is 
overly conservative. 

Other authors have made indire ct measurements of the 
LyC e scape fraction at z ~ 3. iFaucher-Giguere et al.l 
(2008) predict (f^ c c ) ~ 5x 10~ 3 to account for the ioniz- 
ing background inferred from the photoionization rate in 
the Lya forest. However, their estimate is averaged over 
all galaxies down to Ljjv = 0, assumes an intrinsic UV 
to LyC flux density ratio of unity, and is computed using 
an ionizing spectral slope of a s = +0.5; all of these as- 
sumptions decrease the e scape fraction need ed to match 
the photoionization rate. iChen et ahl (|2007l ) use spectra 
of long-duration 7-ray bursts to determine the column 



densities of neutral hydrogen along the sightline though 
the ISM. Their results suggest (f^jf) = 0.02 in the sub- 
L* galaxies that are expected to dominate their sample. 
It is unclear if the difference in the various estimates of 
the escape fraction are due to subtleties in the methods 
or differences in the samples. For example, it could be 
that {fes,p) depends on environment or galaxy type. 

Our sample of NB3640-detected LBGs is small, and 
likely ~ 50% contaminated by non-detections. Thus, 
any differences in the average rest-frame non-ionizing 
UV properties between galaxies with and without leak- 
ing LyC flux may not be evident in our data unless the 
trends are particularly strong. Nonetheless, we inves- 
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tigated the UV properties of the LBGs in our sample 
by cr eating stacks of LRIS spectra (Sha plev et al.ll2003l 
2006) for galaxies with and without NB3640 detections, 
which we show in Figure 

EE We find no significant 
differences between the two samples in either their av- 
erage spectral slopes or interstellar absorption lines, al- 
though the strength of the Lya emission line is on aver- 
age weaker in our LBGs with NB3640 detections. The 
LBGs undetected in NB3640 are also slightly fainter on 
average, consistent with our findings from the stacked 
i?-band photometry (Table [6]). We also estimate the av- 
erage UV slopes of the two subsamples photometrically, 
both through aperture photometry on stacked images in 
the V- and i?-bands, and through averages of the indi- 
vidual fluxes after correcting the V-band magnitudes for 
Lya emission and IGM absorption. With neither method 
do we find significant differences in the average V — R 
colors between the two subsamples. We do note that the 
26 LBGs in our sample are on average slightly bluer in 
G — R (corrected for Lya emission and IGM absorption) 
th an LBGs in the larg e spectroscopic sample presented 
in iShaplev et al.l (|200 3), which highlights the need for a 
larger "average" sample. 

5.1.2. LAEs 

We detected 27 of the 110 main sample LAEs in 
the NB3640 image. Our Monte Carlo simulations 
suggest that 18 - 24 of the detections are uncon- 
taminated, leading to a corrected LAE detection rate 
of 16 - 22%. After applying contamination- and 
IGM absorption-corrections we determine an ensem- 

ble {FuvjFLyc)^ = 2.2±g;§. Most of our LAEs 
with individual detections have flux density ratios in- 
consistent with predictions of stellar population mod- 
els even for young, low-metallicity populations and ne- 
glecting the IGM absorption correction. In contrast, 
the lower limit on the average UV to LyC flux den- 
sity ratio for LAEs without individual NB3640 detec- 
tions ((Fuv I FL y c)o^ s E ' non ~ det > 12.2) is more than a 
factor of 9 higher than the average observed ratio for 

those with detections ((Fuv / FL V c)^ s E ' d&i — 1-3), im- 
plying a dichotomy in the far-UV properties of z ~ 3 
LAEs. 

The implied relative escape fractions flf c > 1 for our 
LAEs detected in NB3640 are difficult to explain. Ex- 
otic stellar population models with top-heavy IMFs and 
extremely low metallicities could in principle reproduc e 
the observed flux density ratios (e.g.. lInoue et al.ll2011h . 
though we caution that there is little additional evidence 
to suggest such models are appropriate for z ~ 3 LAEs. 
As with the LBGs, we find no significant differences be- 
tween the average non-ionizing UV slopes of LAEs with 
and without NB3640 detections, determined either via 
aperture photometry on the R- and V-band stacks, or via 
Lya-corrected averages of the individual V and R mag- 
nitudes. In Figure 1X51 we show the Lya rest-frame equiv- 
alent widths for our LAE sample, estimated from their 
BV— NB3980 colors, as a function of UV continuum (i.e., 

7 We note that the offset of the NB3640 detection from the R 
band position in some LBGs is significant compared to the 1'.'2 slits 
through which the LRIS spectra were obtained. Thus, the stellar 
populations associated with the NB3640 emitting regions may not 
contribute to the observed spectra in these cases. 



i?-band) magnitude. The largest equivalent widths tend 
to be confined to fainter continuum sources, consistent 
with past findings (e.g., lAndo et "all 120061 : iStark et "all 
12010( 1. In order to be detected in NB3640, the faintest 
sources require lower non-ionizing UV to LyC flux den- 
sity ratios (e.g., see Figure ITT]) . Thus, any comparison 
of the distributions of equivalent widths between sources 
with and without NB3640 detections should be confined 
to a relatively narrow range of continuum magnitude. 
The inset in Figure [15] shows such distributions for the 
ranges R < 26.1 and 26.1 < R < 27, demonstrating that 
the median Lya equivalent widths are smaller for LAEs 
detected in NB3640 compared to those without detec- 
tions. 

5.1.3. Faint LAEs 

In addition to our main sample LAEs, we have identi- 
fied 20 "faint" LAE candidates with 26 <NB4980< 26.5 
and BV-NB4980> 1.2. These faint LAEs have larger 
photometric uncertainties (e.g., see Table [2]). As they 
have different selection criteria from the main sample and 
are detected at lower significance, we have not included 
them in our statistical analyses. However, for complete- 
ness we present results for the average properties of this 
sample from our stacking analysis in Table [6] and for in- 
dividual objects in Table fTUl We detect only one object 
in the faint sample in NB3640, while in the main sample 
we detect 14% of the 65 LAEs having W-NB4980> 1.2 
in the NB3640 image. However, the faint LAEs are also 
on average fainter in R and therefore have smaller detec- 
tion limits in Fuv I 'Fl v c- Thus, the data do not allow 
us to address any statistical differences between the two 
populations. 

5.2. Other considerations 

5.2.1. Are we missing low (Fuv / Fl v c) obs LBGs? 

Of the 110 LAEs in our sample, 24% have 
NB3640-E < 1.5, compared to 8% of our LBGs. 
The Lyman break technique specifically selects against 
galaxies with nota bly blue broadband U n — R colors 
(jSteidel et al.l I2003T ) . and while it is possible for a z > 
3.06 galaxy to have significantly higher flux density 
through the narrow NB3640 filter than through the rel- 
atively broad U n -ba,nd, we note that the effective wave- 
lengths of the two filters are very similar. The LAEs 
discussed here, on the other hand, are selected by their 
Lya emission line strength independently of their rest- 
frame UV colors and corresponding (Fuv / Fl v c) values. 
Therefore we must consider whether the very applica- 
tion of the Lyman break technique selects against a sig- 
nificant number of non-LAE galaxies with R magnitudes 
similar to our LBG sample but NB3640— R colors similar 
to those observed in our LAE sample. If such a popula- 
tion exists at z ~ 3, it would appea r as a high-redshift 
tail in samples of BX/BM galaxies (jSteidel et al.l 120041: 
iReddv et al.ll2008h selected to target galaxies with bluer 
U n — R colors at 1.4 < z < 2.5. However, the fraction and 
associated surface density of such galaxies with z > 3.06 
is ex ceedingly small (see, e.g., Figure 1 of IReddv et al.l 
2008), and the incidence of such systems is likely con- 
sistent with foreground contamination. Thus, while the 
nature of the high-z BX/BM galaxies should be explored 
further, it does not appear that we are missing signifi- 
cant numbers of rest-frame UV bright galaxies having 
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Fig. 14. — Comparison of composite rest-frame UV spectra for the six LBGs with LyC detections in the NB3640 filter (black), 
and the 20 objects without detections (red). The composite spectra have been scaled to a common value over 1400 - 1500A. 
In terms of overall spectral shape and strength of interstellar absorption features, the two spectra are very similar. The only 
notable difference is observed in the strength of Lya emission (indicated in greater detail in the inset panel). The spectrum of 
NB3640 non-detections has a Lya equivalent width ~ 2 — 3 times larger than that of the detections. At the same time, the 
centroid of the Lya profile in the spectrum of the detections is shifted towards longer wavelengths. 



NB3640— R ~ due to our selection of z ~ 3 galaxies 
using the Lyman break technique. 

5.2.2. Viewing angle effects 

As with the bulk of our z ~ 3 sample, searches for 
emerging LyC photons from galaxi es at lower redshifts 
have largely yielded null results (iMalkan et al.l 120031 : 
ISiana et all [20071 120101 IBridee et all 12010ft . The small 
fraction of star-forming galaxies with detected LyC emis- 
sion implies that, for ionizing radiation to escape, there 
must be pathways through the neutral ISM from star- 
forming regions that are unusually clear of dust and 
gas. Processes capable of removing or destroying neu- 
tral gas and dust (e.g., feedback, tidal stripping, etc.) 
in actively star-forming galaxies are unlikely to be suf- 
ficiently effective to allow the transmission of ioniz- 
ing radiation over all solid angle. Large spatial vari- 
ations in f^vp are also predicted by simulations (e.g . , 
iRazoumov fc Sommer-Larsenl[2007l : iGnedin et al.ll2008h . 
If LyC flux only escapes over a small fraction of the total 
solid angle relative to the non-ionizing flux, this would 
naturally explain the large spread in the observed far-UV 
properties in our samples. 

We consider a simple model in which most of the UV 
and LyC flux emitted in directions of cleared ISM escapes 
while some UV and no LyC flux escapes in other direc- 
tions. If the projected area of the cleared ISM is compa- 
rable to the size of the LyC-emitting region, the relative 
escape fraction from that region would either be negli- 
gible or near unity, depending on whether the viewing 
angle samples a direction through the cleared ISM. An 
ensemble of galaxies having a single, compact starburst 
would then be expected to exhibit a bimodal observed 
distribution of relative escape fractions. In contrast, a 
single projected view of a galaxy having multiple or ex- 



tended regions of star formation could sample regions 
both cleared of and obscured by neutral ISM, effectively 
averaging out the bimodality inherent in compact sys- 
tems. A wider range of observed values for the relative 
escape fraction would then be possible for such a galaxy. 
Based on the ACS-F814W images in Figures |3] - the 
LAEs in our sample appear to be on average more com- 
pact than the LBGs and, consistent with this simple pic- 
ture, more bi-modal in their LyC properties relative to 
the LBGs. 

If the escape of LyC photons is highly anisotropic, then 
the observed UV to LyC flux density ratio in an individ- 
ual source is unlikely to be representative of the ratio 
of the integrated (i.e., over all solid angle) UV to LyC 
flux densities escaping from that source. For a sample 
large enough to contain a significant number of LyC de- 
tections, however, the ensemble average ratios should re- 
flect the average integrated flux density ratios. Addition- 
ally, under the assumption that all galaxies in each sub- 
sample have similar properties, the detection rate should 
indicate the average fraction of solid angle over which 
LyC photons can escape a typical star-forming galaxy at 
z ~ 3. Of course, if there are evolutionary effects within 
our LBG and LAE samples, the assumption of equiva- 
lence within each galaxy sample is not valid. Careful 
consideration of the stellar populations of the sources is 
required to test this assumption (Steidel et al., in prep). 
For now we note that large differences are not apparent 
in the average UV properties of LBGs or LAEs with and 
without NB3640 detections. 

5.2.3. IGM coherence effects 

When determining the uncertainty in the IGM absorp- 
tion corrections for the ensemble samples, we assume the 
IGM properties along the sightlines to each source are 
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Fig. 15. — LAE Lya rest-frame equivalent widths determined from BV~ NB4980 colors. Sources detected in NB3640 are 
indicated by circled points and open arrows (plotted slightly larger than solid arrows for clarity). Consistent with previous 
studies, we find that the highest equivalent widths tend to be associated with sources that are fainter in the rest-frame UV 
continuum. The inset shows the distribution of equivalent widths for brighter (left) and fainter (right) LAEs; the median 7?-band 
magnitude for sources with i?-band detections is R = 26.1. Sources with and without NB3640 detections are represented by the 
black hashed and red histograms, respectively. 



uncorrelated. If this is not an adequate assumption, the 
uncertainties in our computed escape fractions and IGM- 
corrected emissivities may be underestimated. Indeed, 
studies have found coherence in the Lyct fores t over Mpc 
scales ([Becker et all 12004 ICasev et all 120081) at z ~ 3, 
though correlations are weaker in the higher column den- 
sity systems relevant to our corrections and a detailed 
description of the sizes and co rrelation scale s of forest 
clouds is still lacking (see, e.g., Meiksin 2009). Further- 
more, there may be systematic environmental effects that 
we are not considering. For example, the fact that our 
targets lie in or be hind a large-scale galaxy overdensity 
(jSteidel et al.l feOOO) may cause the IGM to be on average 
more opaque than expected with a higher than average 
incidence of LLSs. Alternatively, a proximity effect due 
to the relatively high density of star-forming galaxies and 
QSOs in the protocluster could in principle diminish the 
IGM opacity on Mpc scales making the sightlines to our 
targets less opaque to ionizing radiation. It should be 
noted, however, that the UV to LyC flux density ratios in 
our LAEs with detections are uncomfortably small even 
before the IGM corrections are made. 

5.3. Implications 

Our results suggest the following picture. Neutral gas 
and dust in the ISM of typical star-forming galaxies at 



z ~ 3, as traced by LBGs, reprocesses most of the non- 
ionizing UV flux, and effectively all of the LyC flux, over 
most solid angle. Over some significant fraction of solid 
angle (~T0 - 20%, as suggested by our corrected de- 
tection rates), however, the neutral gas and dust have 
been cleared or destroyed, allowing both ionizing and 
non-ionizing flux to escape into the IGM in those di- 
rections. Such regions are likely to be smaller than the 
projected UV sizes of typical bright galaxies but may 
be comparable to the sizes of LyC-emitting regions in 
compact sources. In this picture, when detected in LyC, 
compact sources should on average exhibit larger rela- 
tive escape fractions compared to LyC-detected galaxies 
with extended star-forming regions. Sightlines to rela- 
tively faint, low metallicity star-forming sources, such as 
may be selected by our LAE sample, could have lower 
intrinsic UV to LyO ratios further boosting the observed 
non-ionizing to ionizing UV flux density ratios. 

We are left with two difficulties in this picture, how- 
ever. First, the observed LyC flux density relative to 
the non-ionizing UV continuum for the majority of our 
sources with putative LyC detections exceeds predicted 
unattenuated ratios for current reasonable stellar pop- 
ulation models. Second, even if such ratios were pos- 
sible to explain, the implied contribution to the in- 
tergalactic neutral hydrogen photoionization rate from 
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LAEs is in excess of that determined from measurements 
of the opacity in the Lya forest (Bo lton et al.l 120051 : 
iFaucher-Giguere et al.| [2008'). A satisfactory model that 
can reproduce all of the observational evid ence is still 
want ing and follow up spectroscopy (e.g., Ilnoue et al.l 
120111) is needed to determine the nature of the NB3640 
detections associated with our LAEs. 

The LyC photons that escape from galaxies prior to the 
end of reionization cannot be directly observe d, as they 
are absorbed during the reio nization process (|Fan et al.l 
120061 : IRobertson et al.|[2010[ ). Even after the epoch of 
reionization, at redshifts z > 4, the opacity of the IGM 
to ionizing photons due to overdensities retaining sig- 
nificant neutral fractions is sufficient to make direct de- 
tection of escaping LyC flux improbable (|Inoue fc Iwatal 
120081: IVanzella et all I2010al) . Thus, determinations of 
feJc C from direct measurements at z « 3 are likely to 
be the best constraint on the LyC escape fractio n for the 
galaxi es that presumably reionized the Universe. iGnedinl 
(200f) determined that, if the LyC escape fraction for 
low mass (Mhaio < IO^Mq) galaxies is negligible, there 
would be insufficient LyC photons for reionization. On 
the contrary our results for LAEs (while puzzling) im- 
ply that relatively faint galaxies on average have LyC 
escape fractions larger than those of bright galaxies. Re- 
cent data obtained with the Wide Field Camera 3 on 
HST have allowe d for the study of galaxies out to red - 
shifts z ~ 7 (e.g iBunker et alj[20Tot besch et al1[2010h . 
IBouwens et al.l (|2010| ) find a steep faint-end slope to the 
z ~ 7 LF, a = —1.94. Parameterizing the observed evo- 
lution of the LF from z ~ 4 to z ~ 8, they extrapolate the 
LF to even higher redshift and claim that the Universe 
can be reionized by z ~ 6 if the average escape frac- 
tion is (feJb) =0.2. This result is in agreement with 
earlier findings bylOuchi et all (l200l . If (f^ y c c ) = 0.6, 
IBouwens et al.l ([20100 claim reionization is complete by 
z ~ 8 and the predicted Thomson optical depth would 
be within la of that measured for the WMAP 7-year 
data set. Since the luminosity density at these redshifts 
is apparently dominated by relatively faint galaxies, the 
LyC escape fraction in fainter sources such as our sample 
of LAEs is of particular importance. Resolving the dis- 
crepancy between the theoretical minimum and observed 
UV to LyC flux density ratios is thus an important un- 
dertaking for future work. 

6. SUMMARY 

We have presented analysis of the deepest near-UV 
image to date of the SSA22a field. This image was 
obtained through a narrowband filter sampling the 
redshifted Lyman-continuum emission from z > 3.06 
sources, including many galaxies in the z — 3.09 ± 0.03 
SSA22a protocluster. The image contains 26 LBGs with 
spectroscopically-determined redshifts z > 3.06 for which 
our NB3640 filter is uncontaminated by emission long- 
ward of the redshifted Lyman limit. We augmented 
these data with both new and archival deep narrow- 
band and broadband images designed to select LAEs at 
3.06 < z < 3.12. From these images we have identi- 
fied 110 LAE candidates. Our main conclusions are as 
follows: 

1. Six of the 26 LBGs have NB3640 detections with 
AB magnitudes in the range 25.2 to 27.0. Five of 



the detections are offset from the i?-band detec- 
tions by < 1". We have used Monte Carlo sim- 
ulations to determine the probability of contam- 
ination by foreground objects and conclude that 
two to four (Ict confidence interval) of the six de- 
tections are uncontaminated. Twenty-seven of the 
110 LAEs are detected in NB3640 with magnitudes 
spanning 24.7 to 27.2, twenty of which are offset 
from the i?-band detection by < 1". Our sim- 
ulations suggest that 18 to 24 of our LAEs with 
NB3640 detections are uncontaminated. 

2. For our sample of LBGs we find a large range in 
NB3640--R colors and thus observed UV to LyC 
flux density ratios. Furthermore, the subsample of 
LAEs having no individual NB3640 detection has a 
lower limit on the average UV to LyC flux density 
ratio that is almost an order of magnitude larger 
than observed for LAEs with NB3640 detections, 
implying a significant dichotomy in the observed 
far-UV properties of LAEs. We interpret this vari- 
ety in the observed LyC properties of z w 3 galaxies 
as arising from the partial clearing of neutral gas 
and dust over a limited solid angle in individual 
systems. In this picture, ensemble measurements 
of fe S y c c (when corrected for foreground contami- 
nants and IGM absorption) will produce accurate 
average values, while measurements in individual 
systems will be highly biased by orientation. 

3. We find no difference between the average non- 
ionizing rest-frame UV properties for LBGs with 
and without NB3640 detections, with the possible 
exception that the Lya emission line may be weaker 
in those with detections. We also find no difference 
in the average UV spectral slopes for LAEs with 
and without NB3640 detections, although we again 
find evidence for weaker Lya emission in LAEs with 
detections. The similar UV properties of the two 
samples supports the idea that orientation effects 
determine the detectability of escaping LyC flux, if 
such affects are less important for the non-ionizing 
UV compared to the LyC. 

4. After applying statistical corrections for foreground 
contaminations and IGM absorption, and assuming 
an intrinsic UV-to-LyC flux density ratio of six and 
a UV escape fraction of ~ 20 - 25%, we determine 
fLyC ~ 0.1 for our sample of LBGs. The LBGs in 
our sample represent the bright end (L > 0.5L*) 
of the rest-frame UV LF at z « 3. If their av- 
erage UV to LyC flux density ratio is representa- 
tive of LBGs down to 0.1L*, their inferred con- 
tribution to the intergalactic hydrogen photoion- 
ization, Thi, rate is consistent (within the large 
uncertainties) with that measured from the Lya 
fores t opacity ([Bolton et al.l 120051 : iKirkman et al.l 
I200a IFaucher-Giguere et alil2008D . 

5. Half of our LBGs and all but one of our LAEs with 
NB3640 detections have IGM-corrected UV to LyC 
flux density ratios significantly smaller than the 
minimum predicted from stellar population mod- 
els. Our LBGs are drawn from the bright end of 
the LBG LF, L > 0.5L*. If the average LAE UV 
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to LyC flux density ratio that we measure is repre- 
sentative of LBGs over the magnitude range of our 
LAEs, 0.1L* < L < 0.5L*, then the inferred contri- 
bution to the global hydrogen photoionization rate 
Th i exceeds the values measured from the Lya for- 
est opacity at z w 3. Further study is needed to 
determine the nature of the LyC properties of these 
sources. 



detections, such data are needed to confirm the redshifts 
of the LAEs. Finally we note that while our samples of 
LBGs and LAEs are large enough to determine robust 
average properties, it may be that the SSA22a proto- 
cluster is itself unique. Thus we have ongoing projects 
to supplement this work with similar studies in other 
fields. 



Recent measurements of the LBG luminosity func- 
tion at redshifts z > 6 suggest a large value of f^ s v c c 
is needed with substantial contribution from relatively 
faint galaxies in order for the Un i verse to be reion- 
ized bv z ~ 6 (lOuchi et al.1 [200l lOesch et al.1 [200l 
120101 : iBouwens et all l2010f ). While the extremely blue 
NB3640— R colors of our faint LAE sample are difficult to 
explain within current theoretical models for the UV and 
LyC emission properties of star-forming regions, recon- 
ciling the observations with models is an important step 
towards understanding the processes involved in reion- 
ization. To this end, we are augmenting our data sets 
with additional HST imaging and ground-based spectra 
of both LBGs and LAEs with NB3640 detections. In 
addition to providing independent confirmation of LyC 
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APPENDIX 

LBGS AND LAES WITH NO NB3640 DETECTION 

In this appendix we present coordinates, i?-band magnitudes and limits on the observed UV to LyC flux density ratios 
for the LBGs and LAEs that were undetected in the NB3640 image, as well as emission and interstellar absorption 
redshifts for the LBGs. Table[9jlist the LBGs and Table [ID] lists the LAEs. The color W-NB4980= 3.49 corresponds 
to infinite Lya equivalent widths; values near this limit reflect the large photometric uncertainties. 



TABLE 9 

Photometry for LBGs without NB3640 detections. 



ID 


RA 
(J2000) 


Dec 
(J2000) 


7 a 


b 

aba 


R 


F UV c 


C32 


22: 


17:25.63 


0: 


16:12.9 


3 


.301 


3.290 


23.64 


> 29.2 


C30 


22: 


17:19.29 


0: 


15:45.0 


3 


.104 


3.097 


23.70 


> 27.5 


C47 


22: 


17:20.24 


0: 


17:32.5 


3 


.075 


3.065 


23.78 


> 25.7 


D3 


22: 


17:32.40 


0: 


11:33.6 


3 


.086 


3.077 


23.92 


> 22.4 


C35 


22: 


17:20.23 


0: 


16:52.5 


3 


.103 


3.098 


24.06 


> 19.7 


C24 


22: 


17:18.94 


0: 


14:45.4 


3 


.102 


3.091 


24.19 


> 17.6 


Cll 


22: 


17:25.68 


0: 


12:35.4 


3 


.109 


3.099 


24.21 


> 17.3 


C12 


22: 


17:35.29 


0: 


12:47.9 


3 


.118 


3.106 


24.22 


> 17.0 


C4 


22: 


17:38.91 


0: 


11:02.0 


3 


.076 




24.28 


> 16.2 


MD23 


22: 


17:28.01 


0: 


14:29.6 


3 


.092 


3.082 


24.34 


> 15.2 


MD14 


22: 


17:37.91 


0: 


13:43.9 






3.094 


24.49 


> 13.4 


M10 


22: 


17:26.80 


0: 


12:21.3 


3 


.103 


3.095 


24.50 


> 13.1 


C48 


22: 


17:18.58 


0: 


18:16.7 


3 


.090 


3.079 


24.57 


> 12.3 


M28 


22: 


17:31.66 


0: 


16:58.0 


3 


.094 


3.088 


24.75 


> 10.4 


C28 


22: 


17:21.13 


0: 


15:27.7 


3 


.076 




24.87 


> 9.3 


C50 


22: 


17:37.68 


0: 


18:21.2 






3.086 


25.01 


> 8.2 


C26 


22: 


17:39.54 


0: 


15:15.6 


3 


.178 




25.01 


> 8.2 


C15 


22: 


17:26.13 


0: 


12:55.4 


3 


.094 




25.02 


> 8.2 


C39 


22: 


17:20.99 


0: 


17:09.5 


3 


.076 




25.04 


> 8.0 


M14 


22: 


17:39.05 


0: 


13:30.1 


3 


.091 




25.20 


> 6.9 



Lya emission redshift. 
k Instellar absorption rcdshift. 

C Lower-limit for observed ratio of non-ionizing UV and LyC emission with no correction for IGM absorption. 



NARROWBAND LYMAN-CONTINUUM IMAGING 
TABLE 10 

Photometry for LAEs without NB3640 detections. 



ID RA Dec 4980 BV-NB4980 EW a R p^- 

(J2000) (J2000) (A) 



001 c 


22 


17 


32 


40 


o 


11 


34 


\ 


23 


00 


0.74 


21 


23.92 


> 22.4 


UUZ 


oo 
ZZ 


1 1 


QQ 
OO 


on 

yu 


U 


1 1 
11 


Ul 


Q 
O 


OQ 
ZO 


lo 


1 A 1 

1.41 


fti 
01 


O/l OQ 

Z4.Zo 


\ 1 ft o 

> 10. Z 


UU4 


on 

zz 


1 7 


Ofl 

zo 


Ul 


n 
U 


1 -1 
11 


OA 

OU 


n 
U 


OQ 
Zo 


/in 

4U 


l.UO 


Qft 
OO 


O/l Q/1 
Z4.04 


v 1 c o 
> 10. Z 


UUO 


oo 

ZZ 


1 *7 


OO 


Qft 


n 
U 


10 




A 

4 


OQ 
ZO 


1A 
i 1 


O ftQ 
Z.Oo 


4UU 


Z0.00 


\ A ft 
> 4.0 


UUO 


ZZ 


"1 7 


O 1 

Z4 


oU 


n 
U 


1 1 
11 


1 ft 
10 


Q 
O 


OQ 
ZO 


1 


1 ftO 

i.oy 


on 

yu 


\ 07 

> Z ( 




UU I 


oo 

ZZ 




7 
Z / 


7Q 
icy 


n 
U 


1 7 

1 / 


Qft 

oO 


o 

y 


OQ 
ZO 


QA 


1 Oft 

i.yo 


1 Q1 

lol 


\ 07 

> z/ 




008* 


22 


17 


21 


11 





15 


28 





24 


00 


0.95 


30 


24.87 


> 9.3 


009 s 


22 


17 


28 


29 





12 


12 


3 


24 


02 


1.66 


87 


25.84 


> 3.8 


Ull 


oo 

ZZ 


1 *7 


OO 


SO 


U 


1 o 
1Z 


1 A 
11 


y 


O 1 
z 1 


i ^ 

10 


1 KA 

1.04 


/ O 


Oft (17 

ZO.U ( 


\ Q 1 
> O.l 


012 h 


22 


17 


31 


69 





16 


57 


6 


24 


23 


0.97 


31 


24.75 


> 10.5 


013 


22 


17 


27 


18 





16 


21 


7 


24 


23 


1.91 


122 


25.98 


> 3.4 


014 


22 


17 


19 


25 





14 


50 


9 


24 


27 


1.42 


62 


25.82 


> 3.9 


015 


22 


17 


21 


84 





12 


12 


7 


24 


28 


1.07 


37 


> 27 




017 


22 


17 


25 


10 





17 


16 


8 


24 


43 


2.14 


170 


26.22 


> 2.7 


020 


22 


17 


37 


33 





16 


31 


4 


24 


55 


1.44 


64 


25.45 


> 5.5 


022 


22 


17 


19 


68 





11 


49 


4 


24 


58 


1.80 


105 


26.11 


> 3.0 


023 


22 


17 


31 


73 





16 


06 


9 


24 


61 


0.88 


27 


24.91 


> 9.1 


024 


22 


17 


34 


17 





16 


09 


7 


24 


68 


2.11 


162 


26.73 


> 1.7 


026 


22 


17 


18 


96 





12 


00 


8 


24 


70 


1.63 


83 


26.59 


> 1.9 


027 


22 


17 


24 


94 





17 


17 


3 


24 


74 


1.59 


79 


26.25 


> 2.6 


029' 


22 


17 


31 


49 





12 


55 





24 


77 


0.82 


24 


25.38 


> 5.9 


030 


22 


17 


21 


75 





11 


38 


8 


24 


81 


1.34 


55 


> 27 




031 


22 


17 


33 


63 





17 


15 


1 


24 


83 


1.67 


88 


26.37 


> 2.3 


032 


22 


17 


26 


61 





13 


18 


1 


24 


83 


1.17 


43 


26.61 


> 1.9 


033 


22 


17 


37 


50 





14 


08 


3 


24 


85 


1.58 


78 


26.62 


> 1.9 


035 


22 


17 


27 


03 





13 


13 


2 


24 


87 


0.90 


28 


> 27 




036 


22 


17 


22 


25 





11 


55 


1 


24 


89 


1.85 


113 


> 27 




037 


22 


17 


20 


96 





18 


07 


3 


24 


89 


0.97 


31 


25.61 


> 4.7 


040 


22 


17 


31 


93 





13 


08 


5 


24 


92 


1.96 


131 


> 27 




042 


22 


17 


21 


50 





17 


04 


7 


24 


93 


1.04 


35 


25.50 


> 5.2 


043 


22 


17 


21 


65 





12 


23 


4 


24 


98 


1.30 


52 


26.24 


> 2.7 


044 


22 


17 


36 


41 





12 


51 





24 


99 


1.36 


57 


> 27 




045 


22 


17 


35 


97 





16 


30 


2 


25 


03 


1.45 


65 


> 27 




047 


22 


17 


36 


05 





15 


06 


9 


25 


04 


1.96 


131 


> 27 




049 


22 


17 


39 


29 





16 


10 


5 


25 


06 


2.17 


177 


> 27 




050 


22 


17 


24 


56 





15 


56 


8 


25 


08 


> 2.42 


> 258 


> 27 




052 


22 


17 


36 


84 





13 


17 


2 


25 


13 


1.04 


35 


26.63 


> 1.8 


054 


22 


17 


39 


05 





11 


33 


9 


25 


18 


1.85 


113 


26.25 


> 2.6 


055 


22 


17 


35 


80 





11 


50 





25 


20 


1.81 


107 


26.67 


> 1.8 


056 


22 


17 


22 


42 





17 


20 


7 


25 


21 


1.75 


98 


26.86 


> 1.5 


057 


22 


17 


25 


40 





10 


58 


3 


25 


23 


1.07 


37 


26.84 


> 1.5 


058 


22 


17 


19 


61 





15 


38 


4 


25 


27 


1.44 


64 


> 27 




059 


22 


17 


24 


98 





12 


30 





25 


29 


0.85 


25 


25.31 


> 6.2 


060 


22 


17 


28 


19 





11 


17 


1 


25 


30 


1.16 


42 


26.61 


> 1.9 


061 


22 


17 


34 


10 





15 


40 


2 


25 


31 


> 2.19 


> 182 


> 27 




062 


22 


17 


22 


87 





14 


41 


7 


25 


31 


1.26 


49 


26.53 


> 2.0 


063 


22 


17 


23 


32 





15 


52 


9 


25 


31 


1.63 


83 


26.55 


> 2.0 


065 


22 


17 


28 


15 





14 


36 


4 


25 


38 


> 2.12 


> 165 


26.91 


> 1.4 


066 


22 


17 


20 


86 





15 


11 


8 


25 


41 


0.73 


20 


26.64 


> 1.8 


067 


22 


17 


36 


26 





13 


11 


7 


25 


42 


1.13 


40 


26.40 


> 2.3 


068 


22 


17 


18 


37 





17 


26 


1 


25 


44 


1.28 


51 


> 27 




070 


22 


17 


39 


28 





14 


00 


2 


25 


45 


0.78 


22 


25.98 


> 3.4 


071 


22 


17 


21 


61 





12 


20 


5 


25 


48 


1.23 


47 


26.77 


> 1.6 


072 


22 


17 


31 


24 





17 


32 


1 


25 


48 


1.78 


102 


27.00 


> 1.3 


073 


22 


17 


39 


12 





17 


11 


7 


25 


50 


1.98 


135 


26.35 


> 2.4 


075 


22 


17 


22 


97 





11 


25 


8 


25 


51 


> 1.99 


> 137 


> 27 




076 


22 


17 


20 


67 





15 


13 


2 


25 


54 


1.13 


40 


> 27 




078 


22 


17 


37 


68 





16 


48 


3 


25 


56 


1.21 


46 


25.95 


> 3.5 


079 


22 


17 


34 


68 





11 


10 


5 


25 


56 


2.19 


182 


> 27 




080 


22 


17 


35 


95 





13 


43 


3 


25 


58 


1.33 


54 


> 27 




082 


22 


17 


35 


44 





16 


47 


6 


25 


60 


1.53 


72 


> 27 




085 


22 


17 


30 


86 





14 


38 


2 


25 


65 


1.16 


42 


26.94 


> 1.4 


086 


22 


17 


28 


42 





13 


42 


8 


25 


65 


> 1.85 


> 113 


> 27 




088 


22 


17 


38 


45 





13 


18 


3 


25 


74 


1.25 


48 


> 27 




089 


22 


17 


38 


54 





15 


22 


5 


25 


74 


> 1.76 


> 100 


26.59 


> 1.9 


090 


22 


17 


18 


25 





14 


06 


4 


25 


75 


1.30 


52 


> 27 




091 


22 


17 


36 


14 





15 


40 


7 


25 


78 


> 1.72 


> 94 


> 27 




092 


22 


17 


23 


97 





15 


27 


8 


25 


79 


> 1.71 


> 93 


> 27 




093 


22 


17 


27 


48 





13 


57 


5 


25 


80 


0.78 


22 


26.35 


> 2.4 


094 


22 


17 


39 


14 





17 


00 


6 


25 


82 


> 1.68 


> 89 


> 27 




095 


22 


17 


37 


19 





13 


28 





25 


82 


1.13 


40 


26.17 


> 2.8 


097 


22 


17 


27 


11 





14 


08 


7 


25 


85 


0.98 


32 


> 27 




098 


22 


17 


24 


01 





13 


19 


5 


25 


86 


> 1.64 


> 84 


> 27 
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TABLE 10 

Photometry for LAEs without NB3640 detections. 



099 


22 


1 ""7 0/" 1 A f 

17:oo.4o 





13:00.3 


25.87 


> 1.63 


> 83 


> 27 






100 


22 


1 7.on ci 

1 ( :oU.ul 





18:11.6 


25.89 


1.26 


49 


> 27 






103 


22 


17.1(1 A c\ 

1 ( :ly\4U 





15:26. 1 


25.93 


> 1.57 


> 76 


> 27 






104 


22 


1 1 .Q7 

1 ( .0 I .DO 





12:55.5 


25.94 


> 1.56 


> ID 


> 27 






105 


22 


17:35.46 





12:23.9 


25.94 


0.81 


24 


26.81 


> 1 


6 


106 


22 


17:zz\8o 





17:57.8 


25.96 


1.54 


73 


26.16 


> 2 


9 


107 


22 


1 ; :zU.yo 





14:46.7 


25.97 


0.76 


22 


26.57 


> 2 





108 


22 


1 *7 O /I 70 

17:24.78 





17:40.4 


25.97 


1.02 


34 


26.84 


> 1 


5 


109 


22 


1 "7 00 no 

17:2,3.98 





17:57.8 


25.99 


0.77 


22 


26.77 


> 1 


6 


110 


22 


17:19.50 





15:57.6 


25.99 


> 1.51 


> 70 


> 27 






111 


22 


17:31.14 





16:42.9 


26.07 


-s. 1 A A 

> 1.44 


> 64 


> 27 






1 1 o 

11Z 


ZZ 


1 f .oZ. ( Z 


U 


1D.04.Z 


ZO.U / 


\ 1 A A 
J> 1.44 


> 04 


> 2 I 






113 


22 


1 ; :Z4.8U 





13:26.9 


OH (17 

26.0/ 


-s. 1 A A 

> 1.44 


> 64 


> 27 






114 


22 


1 ""7 /i tr r\ 

17:34.50 





14:20.0 


26.08 


> 1.44 


> 64 


> 27 






115 


22 


17:33.46 





17:01.2 


26.09 


> 1.43 


> 63 


> 27 






116 


22 


1 ""7 00 nn 

17:z8.U0 





12:14.2 


26.10 


> 1.43 


> 63 


> 27 






117 


22 


1 v. on no 
1 ; :oy.Uo 





12:01.9 


26.11 


> 1.42 


> 62 


> 27 






119 


22 


1 ; :zo.uo 





12:47.8 


26.18 


> 1.38 


> 59 


> 27 






t on 

1ZU 


ZZ 


1 i .ZD. / D 


u 


lu.oy .0 




s> 1.00 


S> DO 


> z 1 






121 


22 


17:26.44 





15:27.5 


26.30 


> 1.33 


> 54 


> 27 






122 


22 


17:38.19 





14:03.7 


26.30 


> 1.32 


> 54 


> 27 






123 


22 


17:35.06 





17:26.0 


26.31 


> 1.32 


> 54 


> 27 






124 


22 


17:22.80 





17:48.7 


26.32 


> 1.32 


> 54 


> 27 






125 


22 


17:38.02 





14:03.6 


26.32 


> 1.31 


> 53 


26.79 


> 1 


6 


126 


22 


17:19.53 





16:48.2 


26.33 


> 1.31 


> 53 


> 27 






127 


22 


17:36.91 





11:27.1 


26.36 


> 1.29 


> 51 


> 27 






128 


22 


17:23.43 





16:07.4 


26.48 


> 1.23 


> 47 


> 27 






129 


22 


17:22.28 





10:57.9 


26.49 


> 1.23 


> 47 


> 27 






130 


22 


17:32.84 





16:48.8 


26.49 


> 1.23 


> 47 


> 27 







a Lya rest equivalent width estimated from BV — NB4980 color. 

k Lower-limit for observed ratio of non-ionizing UV and LyC emission with no correction for IGM absorption. 

c LBG D3, z em = 3.086. 

d LBG C4, z em = 3.076. 

c LBG MD23, z em = 3.092. 

f LBG C28, z cm = 3.076. 

g LBG candidate C9. 

h LBG M28, z em = 3.094. 

1 LBG candidate M13. 



